RATIONAL SOLUTIONS OF THE SCHLESINGER SYSTEM AND 
ISOPRINCIPAL DEFORMATIONS OF RATIONAL 
MATRIX FUNCTIONS I. 



VICTOR KATSNELSON AND DAN VOLOK 



Abstract. In the present paper we discuss the general facts concern- 
ing the Schlesinger system: the r-function, the local factorization of 
solutions of Fuchsian equations and holomorphic deformations. We in- 
troduce the terminology "isoprincipal" for the deformations of Fuchsian 
equations with general (not necessarily non-resonant) matrix coefficients 
corresponding to solutions of the Schlesinger system. Every isoprincipal 
deformation is isomonodromic. The converse is also true in the non- 
resonant case, but not in general. 

In the forthcoming sequel we shall give an explicit description of 
a class of rational solutions of the Schlesinger system, based on the 
techniques developed here, and the realization theory for rational matrix 
functions. 



NOTATIONS. 



• C stands for the complex plane; 

• C stands for the extended complex plane (= the Riemann sphere): 

C = CUoo; 

• C" stands for the n - dimensional complex space; 

• in the coordinate notation, a point t 6 C" is written as t = (t\, . . . , i n ); 

• C" is the set of those points t — (ti,...,t n ) £ C™, whose coordinates 
ti, ...,t n are pairwise different: C™ = C" \ ( Ui<i,j<n,i#j{* : *» = 

• 9Jtfc stands for the set of all k x k matrices with complex entries; 

• [., .] denotes the commutator: for A, B e M kl [A, B] = AB - BA; 

• I stands for the unity matrix of the appropriate dimension; 

• Ti-idJlk) stands for the set of all rational 9Jt/c-valued functions R such that 
detR(z) £ 0; 

• 'P(-R) stands for the set of all poles of the function R, Af(R) stands for the 
set of all poles of the function i2 -1 ; V(R) is said to be the pole set of the 
function R, J\f(R) is said to be the zero set of the function R. 
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1. THE SCHLESINGER SYSTEM 
The Schlesinger system is a Pfaffian system of differential equations of the form: 



f dQj _ [Qj,Qj 
dU 



(1.1) 



t% tj 



dQ l 
dU 



E 



[Qi j Qj 



l<j<n 1 3 



1 < i, J <n, j, (HTl a ) 
1 < i < n, HHJb) 



where Qi, ■ • • , Q n are square k x k matrix functions of t = (t\, . . . , t n ), t G C™. 
The system of equations can be rewritten in a compact way as 



(1.2) 



vQv = ^2[Q t ,,Qv]vlog(t fi -t„), u = l,. 



where c) is the exterior differential. The Schlesinger system l|l.l(l is over-determined. 
The compatibility (the integrability) condition for this system is of the form 



(1.3) 



A E [Q„,Qu\^iog(t„-t u )) = o, v = i 



The condition 1)1 ..''>H must be fulfilled if Q% (t), . . . , Q n (t) satisfy the equations l|1.2|l . 
This condition expresses the equality V^Qv) — 0, which must be satisfied since 

Let us present the integrability condition (|1.3|l in a more explicit form. Since 
OTlog(t M -t„) = 0, 



(1.4) [Q^Qv^log^-tv)) = 



where we have to substitute the expression of the form (|1.2[) for DQ^ and UQ„ in the 
right hand side of l|1.4fl . Performing the substitutions, we obtain the expressions 



n n 



( L5 ) E E [QajQm].^]^^^-*/.) At)log(t M -^) 



Since lu A lu = for every differential form to, 0log(t A1 — i„) A t)log(i M — i„) = 
0. Therefore, the integrability condition for the Schlesinger system (|1.2H can be 
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presented in the form 



(1.6) [[Q\,Qi*],Qv]*toE(t\-tn)M>tos(tp-tv) 



A=l,/x=l; 



[Q^[Qx,Qu]]Vlog{t x -U)AJ)log(t IJ ,-t v ) = 0, u = l, 



,n. 



At^/j, X^v^^v 



Lemma 1.1. For every holomorphic square matrix functions Q±,..., Q n , the in- 
tegrability condition p. 6(1 holds. 

The proof of this lemma is based on two identities. The first of them is the Jakobi 
identity 

(1.7) [[A,B],C] + [[B,C],A] + [[C,A],B}=0, 

which holds for arbitrary square matrices A,B,C (of the same dimension). The 
second identity is an identity related to the differential holomorphic 1-form uJx.p,'- 

(1.8) UA,„= d l X ~f" , l<A, M <n,A^ At . 

t\ — tfj, 

Since 

(1-9) cj A:M = ^,a, 1 < A, (i < n, A ^ fi, 



n(n_ 1) JS< 
2 



there are different forms uj\ tfl . The second identity is: 



(1.10) uj\^ A w^u + w M ,„ A u v ,x + Wi/,a A uj\^ = 0, 1 < A, fi, v < n, X ^ fj, ^ v. 

To prove this identity we observe that the left-hand side of (|1.1UI) is the sum of 
the expressions that are obtained from lu\^ Alj^^ by all cyclic permutations of the 
indices A, /i, v. Computing the exterior product u>x,/j. A u)^ v : 

dt\ A dta + dtn A dt v + dt v A dt\ 
(L11) ^ U »" = (ti-t^-tj ' 

we see that it is represented as a fraction whose numerator is invariant with respect 
to cyclic permutations of the indices A, fi, v. Therefore, the right-hand side of (|1.10|) 
is equal to 

1 1 1 



(tx-t^it^-Q (tp-tuXtu-tx) {t v -tx){tx-t^) 

x (dtx A dtfj, + dtfj, A dt v + db v A dt\) 

Evidently, the expression in the curly brackets vanishes. The identity (|1.1UI) is 
proved. 

Remark 1.2. The forms u>\ „, 1 < A, a < n, A^«, generate the cohomology 

2ni 

ring if*(C",Z) of the space C™. This result was established by V. Arnold, |Arn| . 
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PROOF of LEMMA O Let us fix v G 1, , . . . , n. First, we transform the first 
summand in (|f .6(1 . Applying the Jakobi identity to the factor [[Qa> Qh] > Qui m the 
first sum, we rewrite (|1.6|) as 

n 

( L12 ) J2 [ ^] > Qx] w a,m a 

n 

+ [[q»,Q\],Qh] 

X=l 5) Lt=l; 

n 

+ X! [[Qa^^jQm]^,!/ A lo^ v = 

A=l,/i=l; 

Since the summation runs over all pairs X,fi 7 we can interchange A and [i in the 
second and third sums of Ijl.f 2|) . Using also 1(1.9(1 and the anti-commutativity of 
the exterior product (uj A a = —a A uj for every differential 1-forms lu and er), we 
rewrite l(f ,f 2(1 as 

n 

(1.13) J2 [ Q»] ' Qa]u;a, m A ^ 

X^fi. X^v,ll^v 

Th 

n 

+ X! [[QjUjQ^iQa]^," Aw^a = 

A = l,/j = l; 

The last equality holds because its left-hand side is of the form 

n 

[[Q^jQi/JjQa] {^a :M A^,„ + w„ : a A uj\^ + u)^ v A uj„ t x} , 

A=l,/i=l; 

where the sums in the curly brackets vanish according to the identity ((f ,1U|) . 

Q.E.D. 

From Lemma Tl .1 l and from Frobcnius theorem it follows that the Cauchy problem 
for the Schlesinger system is locally solvable for every initial data. More precisely, 
the following result holds: 

Theorem 1.3. Let t Q G C", Q", . . . , Q° G Ttk- Then there exist a neighborhood 
V(t°,S) of the point t° , 

V(t°,5) = {te C„ : |ti - i°| < 6, . . . , \t n - t°J < 5}, V(i°, S) c C 

and matrix functions Qi(t), . . . , Q n (t), 

Q v (t):V(t°,8)^Wl k ,v = l,...,n, 



RATIONAL SOLUTIONS OF THE SCHLESINGER SYSTEM 



5 



which are holomorphic in V(t ,8), satisfy the Schlesinger system for t G 

V(t°,S), and meet the initial condition 

(1.14) Q v (?) = Ql,v=l,...,n 

at the point t° . Here 5 is a positive number which depends on t° and Q®,v = 
l,...,n. 

Thus Theorem 11.31 is a local existence result. Actually, a much more stronger 
global existence result holds for the Schlesinger system. It turns out that the 
solution Qi(t), . . . , Q n {t) of the above introduced Cauchy problem for this system 
exists globally as meromorphic (multi- valued) functions in C". We do not discuss 
this more general result in detail - this will lead us too far. However, it is worth 
mentioning that every solution of the Schlesinger system can be represented as a 
fraction of holomorphic functions in C™, where some complex valued holomorphic 
function r(.) on C" appears in the denominator. This is the so-called r-function 
for the considered solution. The singularities of this solution are precisely those 
t G C™ where r(t) = 0. The "zeros" of the function r(.) are responsible for the 
singularities of the solution Qi(.), . . . , Q n {-)- 

2. THE TAU-FUNCTION Of THE SCHLESINGER SYSTEM 

Lemma 2.1. Let matrix-functions Qi(t), . . . , Q n (t) satisfy the Schlesinger system 
(|l.l(l . Then the differential 1-form to: 

n 

(2-1) lu^- tr(Q A (t)Qp(t))9log(*A-tp), 

\.fj,=i 

is closed: 

(2.2) duj = 0. 
PROOF. It is clear that 

1 ™ 

too = -J]tt((}Q A (t))Q^))jlog(t A -g 

1 " 

+ 2 Y tr (QAW(fQ Al (i)))flog(iA-^). 

A,M=1 
A#p 

Interchanging the summation indices A and [i in the first sum (and using the rule 
trAB = trBA), we obtain (after renaming A — > v, v — > A) that 

n 

(2.3) 5w= J tr(Qx(t)(DQ v (t)))l)log(t v -tx). 

Substituting the expression 1)1. 2|l for into l|2.3|l . we obtain the following expres- 
sion for Duj: 

n 

(2.4) £)^= J2 ^{Q\[Q^,QA) 

A=l,/i=l,i/=l 
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where is defined by (|1 .S|> . or 

n 

(2.5) da; = ^2 tr (<2 

A=l l( U=l,!/=l 

The traces in the last formula are invariant under a circular permutation of the 
three indices A, \x, ^. So we can replace cj^,,, A u) v ,x by 3 (w M ,„ A u) Vt x + u) v ,x A u>x,fj, + 
Ux^Aui^v). But this quantity vanishes identically according to the identity H1.10JI . 
Thus, 1)J=0. Q.E.D. 



Since a closed form is locally exact, the following theorem is an immediate 
consequence of Lemma 12. II 

Theorem 2.2. Let holomorphic matrix-functions Qi(t), . . . ,Q n (t) satisfy the 
Schlesinger system HI. If) in some simply connected domain T>,T> C C™ . Then there 
exists a holomorphic non-vanishing function r(i) in V such that 

1 " 

(2.6) 91ogr(t) = - ti(Q x (t)Q^t))d\og(t x -t^) 

Such a function r is unique up to a constant factor: if Ti(t) and r 2 (i) are two 
functions satisfying the relation (|2.6|) . then r 2 (t) = CT\{t), where c ^ is a constant. 

Definition 2.3. A function r related to a solution Qi, . . . , Q n of the Schlesinger 
system Ijl.lll by the equality l|2.6|l is said to be the r-function for this solution. 

The exactness of the differential form (|2.1|) . stated in our Lemma f2. II was estab- 
lished by M. Sato, T. Miwa and M. Jimbo in |SMJ2j (see §2.4 there). In |SMJ2j . the 
potential (|2.6|l was introduced as well as the term "r-function" for this potential. 

Remark 2.4. ^From Lemma f2. II it follows that the r-function for a solution of the 
Schlesinger equation is a holomorphic non-vanishing (multi-valued) function for 
those t where the solution is holomorphic. Actually, the r-function is holomorphic 
on the whole C™. It vanishes for those t where the solution is singular. (See the last 
paragraph of the section^ . In a very general setting, the holomorphy of r- functions 
was proved in |Miwj . 

Remark 2.5. The r-function is a very general concept related to non- linear dif- 
ferential equations. This function was first discovered by R. Hirota as a tool for 
generating many-soliton solutions. (References to the works of H. Hirota can be 
found in bibliography of the book [New]). Let us give examples of r- functions for 
the Korteveg - deVries equation 

(2.7) u t - Quu x + u xxx = 0. 
N - soliton solution 

d 2 

(2.8) u(x,t) = -2—^logA(x,t), 



(2.9) A(x,f)=det (6 jk + 



C j C k 

Vj + Vk 



exp[-(r?3 + nl)t - (r)j + r) k )x] 



j,k=i, 



.N 
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N - phase solution 

d 2 

(2.10) u(x,t) = logtf{ax + bt + c), 

where $ is the Riemann theta function; 
rapidly decreasing solution 

(2.11) u(x, t) = -2— log A [Xt00) (t), 

with A[ I OO )(t) the Fredholm determinant of the Gelfand - Levitan equation. In 
these examples the r-function is A(x, t),i?(oi + bt + c), A[ x oc ^(t), respectively. 

In the paper WNTB , r-function appeared as the spin - spin correlation function 
for the two - dimensional Ising model. In JAIM^ . the Fredholm determinant 
det(7 — XK) of the integral operator K with the kernel K (x, x') = Sln j^~f ^ on the 
union of intervals Ui<y< n [a2j-i, a2j]> Oi < • • • < fl2n, is expressed as the r-function 
r(ai, . . . , a2n) of the Schlesinger system constructed from the operator. (See in 
particular section 7 of |JM MS , where the deformation theory of the Fredholm 

integral operator with the kernel Sln j^~f ' was related to the deformation theory of 
linear differential equations). The interrelation between the deformation theory of 
linear differential operators, the Schlesinger system and its r-function, correlation 
function in statistical physics theory, quantum fields and Fredholm determinants 
are covered in the expository survey [SMJlj . The r-function for general matrix 
linear differential equation with rational coefficients was considered in |JMUj . and 
was studied in great detail in j.TMlj . |.TM2| . 

However, the importance of the r-function and its central role in the soliton 
theory was not understood before the works of the research group from Kioto (M. 
Sato, T. Miwa, M. Jimbo, E. Date, M. Kashiwara, Y. Sato). The work of this 
group, surrounding Mikio Sato, reached its peak in the years around 1981. In the 
paper ,SaSa] , the time evolution of a solution of a non- linear equation is interpreted 
as the dynamical motion of a point of the infinite dimensional Grassman manifold. 
In the papers of E. Date, M. Jimbo, M. Kashiwara and T. Miwa, Lie algebras and 
vertex operators were incorporated into the soliton theory. (See the series of seven 
notes DJKM1 , and the paper DJKM2 ). In this setting, the r-function appears 
in a very general way. The paper |Sa| explains some ideas from jSaSaj in more 
detail. The further development of the approach related to the consideration of 
soliton equations as dynamic systems on an infinite dimensional Grassman manifold 
is done in [SeWij . In particular, in section 3 of |SeWi| the geometric meaning 
of the r-function is explained: this function appears as an infinite determinant 
related to the determinant bundle over the Grassman manifold. In MSW], the 
r-function is defined as the generalized cross-ratio of four points of the Grassman 
manifold. Tau function is discussed in Chapters 7 and 8 of the book [Dick . The 
book |MJD| is a presentation of the algebraic analysis of the nonlinear differential 
equations of the KdV type written for students. In particular, it contains a very 
clear presentation of the circle of problem related to r-functions. Chapter 4 of the 
book New is dedicated to the r-function and its relation to the Painleve property 
and to the Backlund transformations. This non-formally written book covers a very 
widespread circle of problems. 
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In the present work we construct a class of rational solutions of the Schlesinger 
system 1)1. For such a solution, we calculate its t- function explicitly. This r- 
function is equal to the determinant of some (finite) matrix which is intimately 
related to this solution. In the considered special case, r-function is a rational 
function of the variables t±, . . . , t n . 



3. THE POTENTIAL FUNCTION FOR A SOLUTION OF THE 
SCHLESINGER SYSTEM. 

From 1)1. li b) it follows that if the matrix functions Qi(i), . . . , Q n (t) satisfy the 
Schlesinger system in a domain Q gC", then 

Hence, there exists a matrix function V(t) in Q such that 
(3.2) Q.( t ) = ^M, teQ, j = l,...,n. 

Remark 3.1. If the domain g is multiconnected, the function V(t) can be multival- 
ued. In this case we can consider V(t) as single-valued function on the universal 
covering cov((/, to) of the domain g, where to is a chosen distinguished point of g. 

Definition 3.2. The matrix function V(t) which satisfies the condition (|3.2I) is 
said to be the potential matrix function for the solution Qi(t), . . . , Q n (t) of the 
Schlesinger system. 

We have proved that the potential matrix function exists for any solution of 
the Schlesinger system. It is clear that the potential function is unique up to an 
arbitrary constant summand. 

Lemma 3.3. Let matrices Q\(t), . . . , Q n (t) satisfy the Schlesinger system 111.1( 1. 
and moreover the condition 



(3.3) Y, Qi® = 

l<i<n 

for every teg. 

Then the function 

(3.4) u(t) d = f J2 teg, 

l<i<n 

is a potential matrix function for the solution Qj. (£),... , Q n (t) ■ 

PROOF. Differentiating the expression on the left hand side of 13.411 . we obtain 

(3.5) =Qj[t)+ ^ ti--gr-. 

■> l<i<n ■> 

Thus, we have to prove that 

dQi(t) _ n 



(3.6) u 



Ki<n 



at 



j 
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Substituting the expression Ijl.lfl for d< gl into the left hand side of (|3.6|l . we 



obtain 

(3.7) a dQl{t) - V u [9i(t) ' Qj " (<)] tj V ^h9M 

, ^ 9tj 1^-^ ti — tj ti tj 

l<i<n J l<t<n, J l<t<n, J 

= Yl [Qi(*),Qi(*)]= E [Oi(*)»Qi(*] = o. 

l<2<n, l<i<n 

Thus, HSU holds. Q.E.D. 

Lemma 3.4. Let Q be a domain in C™, and Zei matrix functions Qi(t), . . . , 
Q n (t) satisfy the Schlesinger system p.ljl /or t E Assume that the equal- 
ity ()3.3|) ZioZcfe at /east at one point t = t$ £ Q. Then this equality holds at 
every point t € Q. 



PROOF. The equality (|3.3jl will be proved if we prove that 
d_ 



J K l<i<n 

/,From it follows that 

dt j dt j ^ dtj 

l<i<n 3 3 l<i<n ■> 



= _ IQMQM + igiiMiW] = o for t € Q 

l<i<n J l<i<n J 

Thus, (EU) holds. Q.E.D. 
The following lemma is an immediate consequence of two previous ones. 

Lemma 3.5. Let the matrix functions Qi(t), . . . ,Q n (t) satisfy the Schle- 
singer system f or every t 6 Q, where Q is a domain in C™. Assume 
that the condition (|3.3|) holds at least for one point t = to € Q. 

Then the function V(t) which is defined by (|3.4|) is a potential matrix 
function for the solution Qx(t), . . . , Q n (t) . 

The expression X)iiQi(t), which appears on the right hand side of the 
equality (|3.4|) . has a useful interpretation in terms of an auxiliary ordinary 
differential equation. Let t = (ti,...,t n ) G C", and Qi, ■ ■ ■ ,Q n be k x 
k matrices. (For the moment, we consider the point t and the matrices 
Qi, ■ ■ ■ ,Q n as fixed rather than variable.) Let us introduce an auxiliary 
complex variable x belonging to the domain 1 C \ {ti , . . . , t n }. Consider the 



We recall that C = C U oo is the extended complex plane. 



10 



V.KATSNELSON AND D. VOLOK 



linear differential equation 

Ki<n 



dx \ x ~Lj 



with respect to a x k matrix function Y (x) of the variable x. The coefficient 
matrix function V" 1 of the equation (13.911 is holomorphic in the domain 
C \ {ti, . . . ,t n }, and if the condition 

(3.10) ^ Q l = 

l<i<n 

holds, then the residue of the coefficient matrix function at the point oo 
vanishes. Thus, under the condition ()3.1Uj) . the Cauchy problem with the 
initial condition 

(3.11) F(oo) = I 

for the differential equation (|3.9j) in the domain C \ {t\, . . . ,t n } is solvable. 
The solution Y{x) of this Cauchy problem is a function holomorphic in the 
disk R < \x\ < oo, where R = max.,- \tA. The function Y(x) admits a 
Laurent expansion in this disk: 

(3.12) Y(x)=I+ — T> R<\x\<oo- 

l<i<oo X 

Lemma 3.6. Under the condition (|3.1Uj) . the Laurent coefficient Y_i of the 
expansion H3.12|) for the solution Y{x) of the Cauchy problem (|3.9j) - (|3. 

can be expressed as 

(3.13) Y- X = - tiQi- 

\<i<n 

PROOF. Under the condition (|3.1Uj) . the Laurent expansion at oo of the 
coefficient matrix function V is of the form 

(3.14) Y = ( Yl UQi) ■ x~ 2 + o{\x\- 2 ) as x^oo. 

l<i<n X 1 M<i<n ' 

Substituting the expansions (|3.12j) and (|3.14|) into the differential equation 
H3,9|) and comparing the Laurent coefficients, we obtain (|3.13|) . 

Q.E.D. 

Theorem 3.7. Let k x k matrix functions Qx(t), . . . , Q n (t) satisfy the 
Schlesinger system for t £ Q, where Q is a domain in C". Let the equality 
(j3.3j) hold at least at one point t = to € Q, or, what is the same 2 , at every 

2 Lemma l3.4l 
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point t E Q. Assume that for each fixed t E Q , k x k matrix function Y(x, t) 
satisfies the linear differential equation 

(3.15) ^^=fj2 ^rP)-Y(x,t), xeC\{ tl ,...,t n }, 

M<i<n % ' 

with respect to x and the initial condition 3 

(3.16) Y(x,t) lx=oc = I. 
Then the coefficient 4 V(t) of the Laurent expansion 

(3.17) Y(x,t) =/-^ + ^4^ + ^4^ + ... , (max|ti|) < |x| < oo, 

X X A x z i 

is the potential matrix function for the solution Qi(t),... ,Q n (t), i.e. the 
equalities ()3.2|) hold. 

PROOF. This theorem is an immediate consequence of Lemmas 13.31 and 13.61 

Q.E.D. 

The differential equation ()3.15j) plays an important role in what follows. 

4. THE SCHLESINGER SYSTEM AS THE COMPATIBILITY 
CONDITION FOR SOME OVER-DETERMINED SYSTEM OF 
LINEAR DIFFERENTIAL EQUATIONS. 

The result stated in this section is an important step in the study of the 
nonlinear Schlesinger system by linear methods. 

Theorem 4.1. Let Q\(t), . . . , Q n (t) be kx k matrix functions, holomorphic 
with respect to t = (ti, . . . , t n ) in some open domain Q, Q C C". 
Assume that the following system of linear differential equations 

l<j<n J 

(4.1b, W— 9 ^ f*" 1 --"' 

dt k x - t k 

for a square matrix function Y is compatible, that is, there exists a k x k 
matrix function Y = Y(x;t), which is holomorphic, satisfies the equations 
(|4.1[) for x E C, t E Q, x ^ ti, . . . ,t n , and, moreover, is not degenerate: 

(4.2) detY(x;t) + 0, for x E C, t E Q, x ^ t x , . . . ,x ^ t n . 

Then the matrix functions Qi(t), . . . ,Q n (t) satisfy the Schlesinger system 
CEEJ) for t E G. 



■^Due to 13.31 . the value Y(x, t)\ x=ao exists. 

^Here we redenote the Laurent coefficient Y_i(i) by V(t). 
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PROOF. Let us exploit the equalities 

— — -— — 1 < k < n 

dtk dx dx dtk ' ~~ ~~ 

Differentiating eauation (|4,la|) with respect to t^, we obtain 

d_dY = ( ^ Qj ) dY ( _ Q fc y 

Substituting the expression for ^ from 1)4. lb|) into the last equality, we get 

(44) ^-^ = - V QjQfc y | Qk Y V -^*-y. 

^Jx-t^x-tk) (x-t k f ^ x - tj - 



Differentiating equation ()4.1b|) with respect to x, we obtain 

d dY Q k y Q k dY 



dx dtk (x — tk) 2 x — tk dx 

dY 
dx 



Substituting the expression for §^ from ()4,la|) into the last equality, we get 



(4 5) d dY _ Q k y QkQj y 

dxdt k {x-tkf i ^ n (x-t j )(x-t k ) 

Thus, equation 1)4.3)1 takes the form 

dQi 

(4.6) y Ml. = y 

Ijtk 

(Since matrix Y is non-degenerate, we can cancel out the factor V.) 

The numerators in (|4.6|) do not depend on x. Thus, (|4.6|) expresses the 
equality of two rational functions. Comparing the residues of both sides of 
1)4.6)1 at x = tj,j ^ k, we obtain Ijl.ll a). Comparing the residues at x = tk, 
we obtain (fl~Tlb\ Q.E.D. 



Remark 4.2. In the proof of Theorem l4.1l we have not exploited the equalities 

d dY d dY 
— — — — = — — — — , 1 < k,j < n. 

dtk dtj dtj dtk 

They lead to (|l.ll a) and provide no additional information on Qi, . . . , Q n . 
Thus, reversing the reasoning in the proof, we can demonstrate that the 
Schlesinger system (|X. 1|> is not only a necessary, but also a sufficient con- 
dition for the solvability of the linear system (|4.1a|) - (|4.1b|) f° r a non- 
degenerate matrix function Y. 
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Remark 4.3. The result, which we formulated as Theorem 14. 11 is contained 
in the paper (Schlj of L. Schlesinger. (See Section II of this paper, especially 
p. 294. See also p. 324 of the book |Sch2| ). This was the paper [Schlj where 
the Schlesinger equations (jl-ljl originally appeared. In [Schlj . the equations 
(II. lj) were derived as the compatibility condition of some over-determined 
system of linear differential equations. In |Garj . R. Gamier has derived 
the Painleve equation Pyi from the compatibility condition of the appro- 
priate over-determined system of linear differential equations. (This over- 
determined system itself appeared in the paper FuR of R. Fuchs, which, 
in its turn, was based on works of his father L. Fuchs published in Sitzungs- 
berichte der Berliner Akademie der Wissenschaften in 1888 - 1898). Since 



the publication of the works [ZaShj and [AKNSj . representing a solution 
of a nonlinear differential equation as the compatibility condition for some 
over-determined system of linear differential equations has become a com- 
mon wisdom. In AKNS , under the strong influence of the paper [ZaShj . M. 
Ablovitz, D. Kaup, A. Newell and H. Segur, have represented the Korteveg de 
Vries equation, the nonlinear Schrodinger equation, the Sine-Gordon equa- 
tion and a number of other non-linear equations of physical significance as 
the compatibility conditions for the appropriate systems of linear differential 
equations. 

5. FUCHSIAN EQUATIONS AND BEHAVIOR OF THEIR 
SOLUTIONS - GLOBAL AND LOCAL 

A Fuchsian equation is a linear differential equation of the form 

dY ( ^ Qj 

l<j<n 



where t±, . . . , t n are pairwise distinct points of the complex plane {t\ i^ti^ 
■ ■ • t n ) , Qj, j = 1, ... ,n are square matrices of the same dimension (say, 
k x k) , x is the complex variable which belongs to the punctured Riemann 
sphere C \ {t%, . . . , t n }. The points t\, . . . , t n are regular singular points for 
the equation (jB.lj) . The residue of the coefficient matrix ^ -^4- of the 

l<fe<n 3 

equation (|5.1|) at the point oo is equal to — Yl Qk- So, if the condition 

Kk<n 



(5-2) 



Kk<n 



is satisfied, the point oo is a regular point for the equation ()5.1j) . In what 
follows we assume that the condition ()5.2j) is satisfied. 

We consider the differential equation 1)5 .lj) as a matricial one. This means 
that we consider solutions Y(x) of this equation which are k x k matrix 
functions. If xq is a regular point for the equation 1)5. 1|) (that is xq ^ 
ti, . . . , t n ), then for every kxk matrix Yq there exists the solution Y(x;xq,Yq) 
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which satisfies the initial condition 



(5.3) 



Y(x;x ,Y )\ 



x=x 



Yq. 



This solution is a holomorphic function in a disk neighborhood of the initial 
point xq. It can be continued analytically along any path which does not 
pass through the singular points ti, . . . ,t n . Thus the solution of the linear 
differential equation (|5.1|) exists globally and is a (multi- valued) holomorphic 
function on the punctured Riemann sphere C \ {t±, . . . ,t n }. For a point 
x £ C \ {ti, . . . ,t n }, the value of the analytic continuation of the solution 
Y(.; xq,Yq) along a path 7(^0, x) C C \ {t±, . . . ,t n }, which starts at the 
point xq and ends at the point x, can depend on this path 7. If j'(xq,x) 
and ^"(xq,x) are two such paths, homotopic in C \ {tx, ■ ■ ■ ,t n }, then the 
values of analytic continuations of the solution Y(.;xq,Yq) along 7' and 
7" coincide. Thus, the value of the analytic continuation of the solution 
Y(.; xq, Yq) to the point x depends on the homotopy class of the path leading 
from xo to x, along which the continuation is made, but not on the path 
itself. In other words, the solution Y(.,xq,Yq) of the differential equation 
1)5.1)1 with the initial condition 1)5.3)1 is holomorphic function on the universal 
covering space cov(C \ {t\, . . . , t n }; xq) of the space C \ {t\, . . . , t n } with the 
distinguished point xq. 

Definition 5.1. A solution Y{x) of a linear differential equation in a con- 
nected domain is said to be fundamental if detF(x) ^ (or, what is equiv- 
alent, if dety(xo) 7^ for some point xo from this domain). 

If Y(x) is a fundamental solution of a linear differential equation, and Z(x) is 
another solution of this equation then Z(x) = Y(x)C, where C is a constant 
(with respect to x) matrix. In what follows we distinguish the solution Y(x) 
of the differential equation (|5.1|) . which satisfies the normalizing (initial) 
condition 



at a distinguished point xo which is regular for the equation (|5.1j) . 

Let 7(xo,xo) be a closed path, which starts at the point xo and ends 
at the same point xq. We call such a path a loop with the distinguished 
point xq. Let Xg(7) be the point of the universal covering space cov(C \ 
{t\, . . . , t n }; Xq), which corresponds to the homotopic class 7 of the loop 
7(xo,xo). We denote by M 7 the value 1"(xq(7); xo, /) of the solution of the 
equation 1)5. 1J) . normalized by the initial condition 1)5.4)1 . at the point Xq(7). 
The correspondence 7 1— ► M 7 is a representation 5 of the fundamental group 
vr(C \ {ti, . . . , t n }; xq) into the general linear group GL(k; C): 

Myy/ = My My/ for every 7', 7" £ 7r(C \ {t\, . . . , t n }; xq) 



As usual, the product J , 'y" is defined by the composed loop, obtained by traversing first the 
loop, generating 7", and then the loop, generating 7/. 



(5.4) 



Y(x ) = I 
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and 



M E = I, where e is the unity of the group 7r(C \ {t\, . . . , t n }; Xo) 



Definition 5.2. The matrix M 7 is said to be the monodromy matrix for the 

equation ()5.1|) . corresponding to the homotopic class 7. The correspondence 
7 1—5- M-y is said to be the monodromy representation of the equation Q5.1JI . 

The group 7r(C\{ii, . . . , t n }; xo) is finitely generated. Let us choose a special 
system of generators of this group. 
Every simple loop (i.e. the loop 
which is homeomorphic to a circle) 
divides the Riemann sphere C into 
two components. Let 7, be a sim- 
ple loop in the space C \ {t 1 , . . . , t n } , 
with the distinguished point xq, 
such that the point tj and the set 
of all other points {t p } p ^j belong to 
the different components and, more- 
over, the loop 7j goes around the 
point tj in the counterclockwise di- 
rection. (See Figure P). We denote 
the homotopic class of this loop by 
7j as well. It is known that the ele- 
ments 71, . . . , 7„ generate the fundamental group 7r(C \ {t\, . . . , t n }; xq) and 

that the equality 71 7„ = e is the only generating relation 6 . 

For what follows, it is crucial to know how a solution of a Fuchsian differ- 
ential equation behaves in a neighborhood of its singular point. 

More generally, let us discuss the behavior of a solution of the differential 
equation of the form 




Figure 1. 



(5.5) 



dY 
dx 



Q 



t 



where Q is a constant k x k matrix, and $ is a k x k matrix function, which 
is holomorphic in the disk Dt tP , centered at the point x = t, of radius p. 
Solutions of the equation ()5.5|) are holomorphic (in general, multi-valued) 
functions in the punctured disk D t<p \ {t}. 

The equation (|5.5|) can be considered as a perturbation of the equation 



(5.6) 



dY 

dx 



The generating relation in this form holds for the loops, arranged as shown in Figure In 

the general case, this relation is of the form ■y pi -y Pn = e where pi,.. . ,p n is a permutation 

of the indices 1, . . . , n which is determined by the arrangement of the loops. 
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The "unperturbed" equation Q5.fi Jl can be solved explicitly. It has a funda- 
mental solution of the form 

(5.7) Y(x) = (x-t) Q . 

(By definition, {x — t)® = g^ 111 ^' - *)). In general, the function (x — t)Q 
considered as a function in the complex plane is multivalued. It can be 
considered as a univalued holomorphic matrix function on the universal 
covering cov(C\{i}) of the complex plane punctured at the point {t} (or, in 
other terminilogy, on the Riemann surface of log(x — t)). The differentiation 
rule 

(5.8) -^(x-tf = -Q-{x-t) Q , x€cov(C\{t}), 
ax x — t v ' 

holds. (The values of (x — in both sides of Q5.8|) are the same). 

The structure of the solution of the original "perturbed" equation (|5.5|) 
depends on the nature of the spectrum of the matrix Q. 

Definition 5.3. A square matrix Q is said to be non-resonant if no two 

eigenvalues of Q differ by a nonzero integer, or in other words, if the spectra 
of the matrices Q + nl and Q are disjoint for every n G Z \ 0. 

Remark 5.4. It may happen that some eigenvalues of a non-resonant ma- 
trix Q coincide (i.e., the eigenvalues of a non-resonant matrix Q may be 
multiple) . 

Proposition 5.5. If Q is a non-resonant (constant with respect to x) ma- 
trix, and <3?(x) is a matrix-function holomorphic in the disk Dt,p, then the 
differential equation (|5.5[) has a fundamental solution Y(x) of the form 

(5.9) Y(x) = H{x){x-t) Q , x G cov(A, P \ {*})> 

defined on the universal covering cov(-Dj iP \ {£}) of the punctured disc Dt,p \ 
{t}, where H(x) is a matrix function, holomorphic and invertible 7 in the 
disk Dt )P - 

The proof of Proposition and even of a more general statement, is given 
in Appendix. 

In the general case, when the eigenvalues of the matrix Q can differ by 
non-zero integers, the situation is more complicated. 

Proposition 5.6. The differential equation ()5.5j) . with a constant matrix Q 
and a matrix function &(x), holomorphic in the disk D tiP , has a fundamental 
solution Y(x) of the form 

(5.10) Y(x) = H(x)(x -t) L (x -t)®, x G cov(A,p \ {*}), 

where H(x) is a matrix function holomorphic and invertible in the disk Dt,p, 
L is a diagonalizable matrix with integer eigenvalues l\, . . . ,lk, and Q is a 



^The inverse matrix function H 1 (x) exists and is holomorphic in Dt tP 
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non-resonant matrix whose eigenvalues Ai,...,A& are related to the eigen- 
values Ai, . . . , Afc 0/ the matrix Q according to 



(5.11) 



Am — Ap Zr>, 



For a certain choice of the integers l\ > 
6e simultaneously reduced to the forms 

L = T- 1 Diag(Zi,...,/ n )T, Q = T 



1 <p < k. 

■ ■ > l n , the matrices L and Q can 



Diag(Ai ! ...,A n )+C7jr, 



where U is an upper-triangular matrix with zero diagonal. 

Proposition 15.61 is applicable to equation (|5.5|) with the general matrix Q. 
In particular, it is applicable to the equation with the non-resonant matrix 
Q. However, if Q is non-resonant, one can choose L = and Q = Q. 

Remark 5.7. It should be mentioned that the matrix function H(x) in (|5.12|) 
(or in ()5.15j0 is not quite unique. The matrix H(t) = H(x) x= t can be an 
arbitrary invertible matrix commuting with the matrix Q. Of course, the 
solution (|5.9|) (or ()5.1U|0 depends on this choice. In the non-resonant case, 
as soon as this choice is made, the local solution of the equation (|5.5|) at the 
point x = t, i.e. the solution of the form ()5.9|) . is determined uniquely (up 
to the choice of arg(x — t) which leads to the non-uniqueness of (x — tp = 
e Qin(x-t)^ j n w j ia ^ follows , we choose 



H(x) 



x=t 



I. 



In the non-resonant case, the latter equality can be considered as a normal- 
izing condition for the local solution (|5.9|) of (|5.5|) at the point x = t. 

However, in the resonant case, the condition H(x)\ x=t = I do not deter- 
mine the solution of the form (|5.1U|) of the equation (|5.5|) uniquely. The 
following example illustrates this non- uniqueness. Consider the differential 
equation 



where 



dY(x) 
dx 



Q 



9- . Y(x) 



This equation, which is of the equation (|5.5j) with t = 0, <J>(x) = 0, has the 
solution 

~x 0" 



Y x {x) = 

which is of the form (|5.1()j) with H(x) 
t = 0. The matrix function 



1 



1 
1 



L 



1 




, Q 








Yo(x) 



or, more explicitely, Y2{x) 



18 



V.KATSNELSON AND D. VOLOK 



is another solution of the same differential equation. However, the function 
Y<z(x) is also representable in the form 



Y 2 (x) 



i.e. in the form (|5,10|) . with H(x) 



1 x~ 




X 


0" 


1 







1 



L 



,Q 



,t = 0. 



Thus, the considered differential equation possesses two different solutions 
Y\{x) and Y2(x), both these solutions are of the form (|5.10|) . with the same 
L and Q, but with different H(x), and both these factors H(x) satisfy the 
condition H{x)\ x= q = I. 

Thus, in the resonant case, the condition H(x)\ x=t = I does not uniquely 
determine the local solution (|5.10)1 of the equation (j5.5j) corresponding to the 
singular point x = t, even sough the matrices L and Q are specified. 

Remark 5.8. There is another essential difference between the non-resonant 
and the "resonant" cases. If the matrix Q is not non-resonant, then the 
matrix Q, which appears in (|5.1U|) . depends on the value = &(x) x= t 
of the "perturbing" (see ()5.5l0 function at the point t. In particular, the 
Jordan structure of the matrix Q may depend on the matrix This 
Jordan structure, together with the eigenvalues of Q, determine the multi- 
valued behavior of the matrix function (x — t)®, and finally, the multi- 
valued behavior of Y(x) in Dt p . Thus, in the "resonant" case the multi- 
valued behavior of the solution Y(x) depends not only on the matrix Q (the 
coefficient of the "leading term") but also on the "perturbing" matrix $. As 
can be seen from (|5.9|) . in the non-resonant case the multi- valued behavior 
of the solution Y(x) depends only on Q, but not on 

Remark 5.9. Of course, for every matrices A and B, the matrix e A 
be presented as e c : 



e B can 



where 



C 



log(e A -e B ), 



where an arbitrary value for the matrix log(e • e ) can be chosen. However, 
the case of exponential functions is different. In general, a matrix function 
e AC . e -BC f £ (with constant matrices A and B) is not representable as e c< ^ 
with a constant matrix C. In particular, in the resonant case the solution 
(|5.1U|) of the differential equation (|5.5|) is not, in general, representable in the 
form Y(x) = H(x) ■ (x — t) R with a holomorphic invertible matrix function 
H and a constant matrix R. 

Propositions 15.51 and l5~6l go back to L. Fuchs and G. Frobenius. The sys- 
tematic study of linear differential equations with analytic coefficients in 
the neighborhood of a singular point was initiated by Lazarus Fuchs in 
1866 and followed by G. Frobenius in 1873. They discussed the "scalar" 
differential equation of the n th order. Matrix differential equations in the 
complex domain were considered later. The statement which we formulated 
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as Proposition 15.51 can be found in a lot of excellent textbooks now. See, for 
example, |CoLe| . Theorem 4.1; |Gant| . Chapter XV, §10; |Hartj . Chapter 
IV, Corollary 11.2; jlKSYj . Chapter I, Theorem 2.3. However, it is difficult 
to point out the source, where this statement first appeared precisely in the 
same form as in Proposition 15.51 The earliest reference we know is |Ra| (see 
Section 3 there). The statement which we formulated as Proposition 15.61 
can be found in [Gant], Chapter XV, §10, in HsSi , Chapter V, Section 5; 
|Sibj . Chapter III, Theorem 3. 9. 4, page 89. The earliest of this references is 
Gant (the first edition of Gantmacher's book appeared in 1954). However, 
Y. Sibuya, |Sibj . page 89, refers to Hukuhara's book |Hukj (Theorem 41.1, 
page 143) published in 1950. Unfortunately, the book |Hukj is unavailable 
to us. 

Let us return to the Fuchsian equation (|5.1|) . We consider this equation 
in a disk neighborhood Dt ,p OI the singular point tj, which does not contain 
other singular points t p ,p ^ j. Let us choose and fix a point Xj of the punc- 
tured disk D tj:P \ {tj}- We also choose and fix a path ctj in C \ {t\, . . . , t n }, 
which starts at the distinguished point xq and ends at the 







Figure 
2. a: 
path ctj. 



Figure 
2. b: 
path [i. 



Figure 
2. c: 
path a. 



point Xj. (See Figure[2Ja.) We need these Xj and ay to distinguish a certain 
analytic continuation Y a . of the solution of 1)5.1(1 normalized by (|5.4[) into 
the punctured disk Dt jtP \ {tj}. Let x be an arbitrary point of the punctured 
disk Dt hP \ {tj}, and fi be an arbitrary path in Dt , p \ {tj}, starting at Xj 
and ending at x. (See Figure EJb.). So, the pair (x, fj) can be considered 
as a point of the universal covering cov(D tjtP \ {tj};xj). We construct the 
composed path a = fi ■ ctj. The path a starts at the point xq and ends at 
the point x. (See Figure[21c.) We continue analytically the solution Y along 
the path a from the point xq to the point x. In this manner we distinguish 
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the continuation Y aj {x) of the solution Y(.), normalized by the condition 
(|5.4|l . from a neighborhood of xq to the punctured disk D tj , p \ {tj}. The 
function Y a . (.) can be considered as a multi-valued function in Dt , P \{tj}, or, 
better to say, as a single- valued function in the universal covering cov(Dt. tP \ 
{tj};xj) with the distinguished point Xj. The notation Y aj reflects how the 
path ctj is involved in the process of the analytic continuation. In fact, this 
process involves two steps. On the first step, the solution is continued along 
the path ctj from a neighborhood of the initial point xq to a neighborhood of 
the point Xj, where the point Xj is considered as a distinguished point of the 
punctured disk D tj:P \ {tj}. On the second step, the solution is continued 
along paths which are contained in Dt jtP \{tj} from the neighborhood of the 
point Xj to points of cov(D tjtP \ {tj};xj). 

The matrix function Y a \.) is a fundamental solution of (|5.1|) . Proposi- 
tions 15 . 5ll5 . 6l provide another fundamental solution of (|5.1(l . As a preliminary 
we assume that the matrix Qj is non-resonant. According to Proposition l5.5l 
there exists the fundamental solution Yj{x) of ()5.1j) . which is representable 
as 

(5.12) Yj(x) = H J {x){x-tj) Q J, x e cov (D tj , p \{tj},Xj), 

where Hj{x) is a matrix function, holomorphic and invertible in the whole 
(non-punctured) disc Dt j)P . Under the normalizing condition 

(5.13) Hj(x) lx=t . = I, 

the matrix function Hj(x) is unique. The matrix function (x — tj)^ j is 
defined up to a right factor of the form exp[2TriQj]. To avoid this non- 
uniqueness, we choose a certain value 

(5.14) Oj = avg( Xj -tj), 

and fix this choice. Under the agreement (|5.14j) . the matrix function (x — 
t)®J is defined uniquely on cov (D tj:P \ {tj}, Xj). Thus, if the matrix Qj 
is non-resonant, the local solution Yj{x) of the differential equation ()5.1|) 
corresponding to the singular point tj can be distinguished uniquely by means 
of the normalizing conditions (|5.13|) and (|5.14|) . 

If the matrix Qj is resonant, the situation is different. According to 
Proposition 15.61 there exists the fundamental solution Yj(x) of ()5.1I) . which 
is representable as 

(5.15) Yj(x) = H j (x)(x-t J ) L >(x-tjf>, x£cov(D tj , p \{tj}), 

where the matrix Hj(x) has the same properties as before, Lj is a diagonal- 
izable matrix with integer eigenvalues (so, the matrix function (x — tj) L i is 
single- valued in Dt, P \ {tj}), and Qj is a non-resonant matrix, whose eigen- 
values differ from eigenvalues of the matrix Qj by some integer numbers. 
However, even though the normalizing conditions (|5.1«-{|) and (|5.14fl are as- 
sumed and the matrices Lj,Qj are chosen, the matrix function Hj{x) in 
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(|5.15j) ) is determined non-uniquely 8 . Thus, if the matrix Qj is resonant, we 
may deal with some local solution Yj(x) of the form (|5.15j) corresponding to 
the singular point tj of the differential equation (|5.1|) , but in general, we can 
not distinguish naturally such a solution. 

Definition 5.10. Assume that the matrix Qj is non-resonant. The solu- 
tion Yj of the Fuchsian equation (|5.1j) representable in the form (|5.12j) near 
{tj}, with Hj satisfying the normalizing condition (|5.13j) . is said to be the 
normalized local solution of (|5.1j) corresponding to the singular point {tj}. 

Let us choose some local solution Yj of the equation (j5.1j) corresponding 
to the singular point tj. If the matrix Qj is non-resonant, it is natural to 
choose the normalized local solution. If the matrix Qj is not non-resonant, 
we choose some local solution of the form (|5.10j) and fix such a choice. Since 
both solutions, Y x ^ a . and Yj, are fundamental, the equality 

(5.16) Y Xj>aj (x) = Y j {x)C, 

holds, where C is a constant (with respect to x) non-degenerate matrix. 

Definition 5.11. The equality (|5.16|) is said to be the connection relation. 
The matrix C from (|5.16j) is said to be the connection coefficient. 

According to Definition I5.11( the connection coefficient C depends on 
the choice of the point Xj € D tjtP \ {tj} and on the choice of the path ay, 
connecting the points xq and xy. C = C Xji0lj . Now we will discuss this 
dependence. 

Let a'j and a" be two paths, both of which start at the point xq and end 
at the point Xj, ot'j,a" C C \ {ti, . . . , t n }. It is clear that a'- = a'j ■ 7, where 
7 is the loop with the distinguished point xq, which is composed of the path 
a'j and of the path (a'j) -1 (the latter is the path a'j which is traversed in 
the opposite direction: from the point Xj to the point xq): 7 = (aj) _1 ■ ct". 
The paths a' and a" , connecting the point xq with the point x, are defined 
as the compositions a' = fj, ■ a'j and a" = fi ■ a'j, where, as before, fi , 
[i C D tjiP \ {tj}, is a path connecting the points Xj and x. Thus, a" = a' ■ 7, 
and 

Y a n(x) = Y a ,.(x)M 1 , 

where M 7 is the monodromy matrix corresponding to the loop 7. On the 
other hand, 

Y a >.(x) = Yj{x)C Xj)a ,., Y a >>(x) = V ; ;./•)( 

Therefore, 

(5.17) C xj,a'! = C Xj ^M y , where 7 = (ay)" 1 • a". 

^See Remark 15.71 In what follows, such an arbitrariness in the choice of the local solution 
does not play any essential role. We need the notions of the local solution and of the connection 
coefficient only to illustrate and motivate the notion of the isoprincipa] deformation which will be 
introduced further, in section 171 
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Let us discuss, how the connection coefficient C depends on the point 
Xj. Given two points x'px'- £ Dt jlP \ {tj} and a path a'j C C \ {t\, . . . , t n }, 
connecting xq and x'-, there exists a path a" C C \ {ti, . . . , t n }, connecting 



xq and x"-, such that 
(5.18) 



Cry. 



Cx'^a" 



Indeed, let v be a path, connecting x'j and x", v C Dt j)P \ {tj}. We define 
the path a" connecting xq with x" as the composition a" 



v ■ a 



The 

pairs (x'-,a'-J and (x'-,o/-) allow us to distinguish the solutions Y a i and Y a /i 

J J J J 3 3 

of the equation (|5.1|) in the punctured neighborhood D tjtP \ {tj}. To obtain 
values of these solutions at a point x £ D tjtP \ {tj}, we have to choose paths 
fi'j and fj,", connecting the points x'- and x'L respectively, with the point x. 
If we make this choice coherently, we obtain the equality Y a i (x) = Y a n (x) 
which implies (|5.18|) . Such a coherent choice can be done in the following 
way. If // is a path in Dt jlP \ {tj}, connecting x'j with x, we define the 
path /i", connecting x" with x, as the composition fj," = fj,' ■ v~ x . To obtain 
the values of the solutions Y x > a / and Y x n a n at the point x, we have to 

j ' 3 j ' 3 

continue analytically the solution Y , which is normalized by (|5.4|) . from the 



fj,' ■ ap- 



point xo to the point x along the paths a' 
respectively. However, these two paths are homotopic. 



and a 



a. 



To see this, we 



• v) ■ a' and that the path 

Y n 



remark that a" = (// • v ) • [y ■ a') = /i' 

u^ 1 ■ v is contractible to the point x\. Hence, Y a i (x) = Y" a »(x). Thus, in 

** 3 3 

some sense the connection coefficient does not depend on the choice of the 
point Xj G D t . >p \ {tj}. 

Let us calculate the monodromy matrix cor- 
responding to the loop "going around" the 
point tj. We relate such a loop 7 Qj to every 

path ctj, otj C C \ {ti, . . . , t n }, which starts at 
the point xo and ends at a point Xj belonging 
to a punctured disk D tj , p \ {tj}. The radius 
p must be chosen so small that the disk does 
not contain any other point t p ,p ^ j. Let [i 
be the circular arc, centered at the point tj, 
of the radius e = \xj — tj\, which starts at the 
point Xj, ends at the point tj + (xj —tj)e 2m and 
is oriented counterclockwise, making precisely 
one turn around its center tj. We construct the 
loop 7 a . as the composition of the path a,-, of 



the arc /j, and of the path a- 
(see Figure I3J). 



7 



[l-OLj 




Definition 5.12. The loop 7^ is said to be the 
loop generated by the path ctj. 



Figure 3. Loop j ar 
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Lemma 5.13. The monodromy matrix M 7a . along the loop "f aj generated 
by a path ctj can be represented as 

(5.19) M 7a . = (Y Qj .(x,))-%(^ + (s; - t 3 )e 2m ), 
where Xj is the endpoint of the path ctj . 

PROOF. For linear differential equations, the following semigroup property 
holds: Let 7 be a path starting at a point a and ending at the point c, and 
let b is an "intermediate" point of the path a, that is the point which is 
located "between" the points a and c. Then Y(c,a) = Y(c,b)Y(b,a), where 
for £,77 G a, Y(rj,£)) is the value Y(x,£)\ x=v of the solution of the equation 
satisfying the initial condition Y(x, £)|a:=£ = ^- Let us take the loop 7^ 
as such a path 7. The point a is the point xq considered as the starting 
point of the loop, the point c is the same point xq, but considered as the 
end point of the loop. The point b is the point tj + (xj — tj)e 2m . From 
the semigroup property it follows that Y(b, c)Y(c, b) = I. It is clear that 
Y(b,a) = Y aj (tj + [xj - tj)e 2 ™), Y{b,c) = Y a .( Xj ), and Y(c,a) = M^. 
Q.E.D. 

From (|5.16j) it follows that in the non-resonant case, 
and 

Yxj,ctj(tj (jEj ) — Hj(xj)e ^{xj tj) ^Cqj-. 

(The function Hj(x) is single- valued in D tj:P , and the function (x — tjP 3 ' 
receives the factor e 27n< ^, when x traverses the arc \x once, in the counter- 
clockwise direction.) Therefore, according to Lemma 15.131 

(5.20) M 7j = C-^Xj - tjY^e^ixj - t^C aj . 
Finally, we obtain 

Lemma 5.14. The monodromy matrix M 7a ., corresponding to the loop ^f a . 
generated by a path ctj can be represented as 

(5.21) M 7j = C~] exp[27riQ,]C aj 
in the non-resonant case, and as 

(5.22) M 7j = C'j exp[27rzg,]C aj 

in the resonant case, where C a . is the connection coefficient, corresponding 
to the path ctj. 

(The reasoning in the resonant case is analogous. One needs to use the 
equality (|5.15|) . taking into account that the matrix- function {x — tj) L i is 
single- valued in the disk Dt 3 ,p)- 
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For what follows, it is very fruitful to rewrite the representations (|5.12j) 
- (jnHSJ), (fS"T5j) - (jSJHl) for a "local" solution V (|5~T|) . as well as the ex- 
pressions (|5,21j) and (|5.22j) for the monodromy matrices, in a different form. 
The respective transformations are trivial, but nevertheless very useful. The 
form in which we present the local solutions hints at how to introduce the 
notion of the isoprincipal deformation of the Fuchsian differential equation. 
This notion will be of crucial importance to us. 

Let D tj<p be a disk neighborhood of the singular point tj of the differential 
equation ()5.1|) . which is small enough, so that the point tj is the only singular 
point which is contained in this neighborhood. In this neighborhood lives a 
"local" solution Yj of the equation ()5.1j) . which in the non-resonant case is of 
the form (|5.12[) . and in the resonant case is of the form (|5.15j) . However, this 
"local" solution Yj does not, in general, satisfy the initial condition (|5.4|h To 
be precise, we should speak about the analytic continuation of this solution 
from the neighborhood Dt^ p to the point xq. To discuss such an analytic 
continuation, we have to choose a point Xj of the punctured neighborhood 
Dt Jt p \ {tj}, and a path ay, ctj C C \ {t\, . . . , t n }, leading from the point 
xq to the point Xj. The solution Y a . which satisfies the initial condition 
(|5.4j) at the distinguished point xo and which is continued analytically to 
the neighborhood Dt jjP along the path aj is of the form 

(5.23) Y aj (x) = Hj(x)(x - t,j) Q >C aj , x £ cov (D tj , p \ {tj}, Xj ) , 
in the non-resonant case, and of the form 

(5.24) Y a .{x) = Hj{x){x-tj) L J(x-tj)QiC aj , x G cov (D tj:P \{tj}, Xj ), 

in the resonant case, where C aj is the connection coefficient corresponding 
to the choice of the homotopic class of the path aj and the choice of the 
arguments &rg(xj — tj), see ()5.14j) . (|5.12j) and (|5.15|) . 
We introduce the matrices 

(5.25) A a . = C a ^QjC aj , H a . = HjC a . 
in the non-resonant case, and the matrices 

(5.26) A a . = C a ^QjC aj , B a . = C a LjC a ., H a . = HjC aj 

in the resonant case. Propositions 15.51 and 15.61 can be "developed" in the 
following way. 

Theorem 5.15. If the matrix Qj is non-resonant, then the solution Y a . of 
the differential equation (|5,1[) that satisfies the initial condition (|5.4|) at the 
distinguished point xq and is continued analytically along a path aj into a 
small punctured neighborhood D t . p \ {tj} of the singular point tj, is repre- 
sentable in the form 

(5.27) Y aj (x) = H aj (x)(x - tj) Aa i , x € cov (D tj , p \ {tj}, Xj ) , 



I.e., for a solution in a neighborhood of the singular point tj. 
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where A a . is a matrix, similar to the matrix Qj, and H aj is a matrix func- 
tion, holomorphic and invertible in the entire neighborhood D t . p (including 
the point tj). 

Theorem 5.16. If the matrix Qj is resonant, then the solution Y a . of the 
differential equation (|5,1|) . that satisfies the initial condition (15,4(1 at the dis- 
tinguished point xq and is continued analytically along a path ctj into a small 
punctured neighborhood Du,p \ {tj} of the singular point tj, is representable 
in the form 

(5.28) Y a .(x) = H aj (x){x-tj) Za i{x-tj) Aa 3 , x G cov (D tj>p \{tj},Xj), 
where Z a . is a diagonalizable matrix with integer eigenvalues A aj 
is a non-resonant matrix, whose eigenvalues X p are related to the eigenvalues 
X p of the matrix Qj by the relation 

(5.29) Xp = X p -lj, l<p<k, 

and H a . (x) is a matrix function, holomorphic and invertible in the entire 
neighborhood (including the point tj). For a certain choice of the integers 
h I'm the matrices Z aj and A a . can be simultaneously reduced to 

the forms 

Z aj = T^Diag (li , . . . , l n ) T a . , A aj = T" 1 (Diag (Ai , . . . , A n ) + Uj T aj , 

where U is an upper-triangular matrix with zero diagonal. 

Remark 5.17. If a'j is another path with the same starting point xq and 
endpoint Xj, then, in the non-resonant case, the matrices A a , and H a from 
(|5.27|) are transformed according to the rule 

(5.30) A a , = M~ l A a . M 7 , H a , (x) = H a . (x)M 7 , 

and in the resonant case, the matrices Z a . , A aj and H a from ()5.28|) are 
transformed according to the rule 

(5.31) Z,j =M~ 1 Z a .M J , A a <. = M~ l A a M 1 , H a ,, [x) = H aj {x)M 1 , 

where M 7 is the monodromy matrix corresponding to the loop 7 = {aj)^ l a'-. 
In particular, if the paths ctj and a'j are homotopic then 

Z a i. = Z a ., A a >. = A a ., H a >.(x) = H aj (x). 

Theorem 5.18. Let M 7q . be the monodromy matrix for the differential 
equation \5.1\) corresponding to the loop 7 Qj . , which goes around the singular 
point tj and is constructed from the path ctj, connecting the initial point xq 
to a small neighborhood of tj: 7^ = (<x,-) _1 ■ fj, - ctj (see Figure^- Then this 
monodromy matrix can be expressed as 

(5.32) M la . =exp(27ra aj ), 

where the matrix A aj is the same as in the exponent in \5.27\) , if Qj is 
non-resonant, and the exponent in \5.28\) . if Qj is resonant. 
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Theorem 5.19. Let M^. be the monodromy matrix for the differential equa- 
tion \5.1\) corresponding to a loop jj with the distinguished point Xo which 
makes one turn counterclockwise around the singular point tj and makes no 
turns around the other singular points t p ,p = 1, . . . ,t n ,p ^ j. Then, what- 
ever the matrix Qj is, non-resonant or not, the spectra spec M 7j and spec Qj 
are related by 

(5.33) spec M 7i = exp {2iri spec Qj}, 

where the latter equality holds "with multiplicities" , that is the multiplicity 
of the eigenvalue \x of the matrix M lj and multiplicities k\ of the eigenvalues 
A of the matrix Qj are related by = ^ k\, were the sum is taken over all 
A such that exp{A} = \i. 

PROOF. This theorem is the immediate consequence of Theorem 15.181 of 
the spectral mapping theorem for matrices (from (|5.32|) it follows that 
specMy^ = exp {2ni spec A a . }), and of Theorems 15.151 or 15.161 which relate 
the spectra of the matrices A aj and Qj. However, here we can do without 
Theorem l5.16l restricting ourselves to non-resonant matrices Qj only. Every 
matrix can be approximated by a non-resonant one. The spectrum spec Qj 
depends on Qj continuously. The monodromy matrix M 7j . depends on Qj 
continuously as well. Q.E.D. 

6. DEFORMATIONS OF FUCHSIAN DIFFERENTIAL 

EQUATIONS 

A deformation of a Fuchsian equation is a family of Fuchsian differential 
equations depending on parameters: 




If the coefficients Q p (a) depend on the parameters a = (ai, . . . , a m ) in some 
open set of the space C m holomorphically, such a deformation is said to be 
a holomorphic deformation of a Fuchsian differential equation. The study 
of linear differential equations that depend on parameters holomorphically 
was started by Lazarus Fuchs, see |FuLlj - |1 uL4| . This study was continued 
by L. Schlesinger, see |Schlj - |Sch3j . and by Richard Fuchs, see |FuR| . In 
particular, L. Schlesinger has considered linear equation, where loci ti, . . . , t n 
of singular points serve as the parameters a. In other words, L. Schlesinger 
has considered deformations of Fuchsian differential equations of the form 




where the the coefficients Q p {t\, . . . , t n ),p = 1, . . . , n, are k x k matrix func- 
tions which are defined and holomorphic for t from a domain T> of the space 
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In the further we assume that the condition 



l<p<ra 



Q p (t) = o, tev, 



is satisfied. 

In what follows we study the normalized by (|5.4|) solution Y(x,t)) of the 
family (|6.1|) of Fuchsian equations as a function of x and t. In particular, 
we consider the monodromy matrices M 7j . for the equation (|6.1|) . as well 
as the matrices A aj and functions H aj (x) from the representations (|5.27|) 
and (|5,28|) for various t\, . . . ,t n : M- Jj = Mj.(t), A a . = A aj (t) and H a . = 
H aj (t), intending to study how these matrices depend on t = (pi, . . . ,t n ). 
Here a problem arises. The matrices M 7j ., A aj , H a . depend on the ho- 
motopy classes of 77 and ctj, which are considered on the punctured Rie- 
mann sphere C\ {t\, . . . , t n }. However, for different t = (t\, . . . , t n ), the sets 
C \ {ti, . . . , t n } are different. Thus we need to explain what does it mean 
that, for the deformation (|6.1|) of the Fuchsian equation, the monodromy 
matrices M 7 ., or the exponents A a . from (|5.27|) (or H5.28j0 are holomorphic 
with respect to t, or what does it mean that these matrices do not depend 
on t. The latter problems are of the local nature. 

Let a point t° = (t®, . . . , G C" be fixed, t® ^ xq, j = 1, . . . , n, where 
xq € C is the distin- 
guished point (the same 
point which appears in the 
initial condition (JH3J, and 
which serves as the dis- 
tinguished point for the 
fundamental group 7r(C \ 

,0; x o) )• Choose 
and fix simple loops = 
l,...,7 n , each of which 
starts and ends at the point 
xq, the loop 7j contains 
the point t® "inside" and 
goes around this point coun- 
terclockwise, and all other 
points {tp} p ^j lie "outside" 
the loop 7j, j = l,...,t n . 
Let 5 > be a small num- 
ber such that the disk D t o s = {z G C : \z — t®\ < 6} does not intersect 

the loop jj,j = l,...,n (and hence, is contained inside the loop 7j). Then 
for every t = (t\, . . . ,t n ) such that tj G D t o s , the loops = 1, . . . ,7„, 

generate the fundamental group 7r(C\{ti, . . . , t n }; xoj). See Figured where 
the disks D t q $ are plotted in the light gray color. 




Figure 4. 
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Assume that the point t° belongs to the domain ~D of the space C™, where 
the matrix functions Qj(t) are holomorphic, and that 5 is so small that the 
polydisk 

n 

(6.2) D(t°,5) = JjD^ where = {t p : \t p - t° p \ < 5}, 

p=l 

is contained in T>: D(t°,8) C T>. Then, for x G jj, the coefficient matrix 
Yl "J~F °f the linear equation ()6.1|) is holomorphic with respect to t in 

D(£°,<5). From the standard results on dependence of solution of a linear 
differential equation on parameters follows 

Lemma 6.1. Under above stated assumptions and notation, the monodromy 
matrices M 7 . (£), j = 1,... ,n, of the holomorphic family (|6.1|) of Fuchsian 
equations are holomorphic with respect to t from the polydisk D(i°,<5). 

Now we investigate in more detail the local solution Yj(x, t) of the equation 
(|6.1|) which corresponds to a neighborhood of the point tj. We have already 
considered such a local solution (see Propositions 15.51 - 15.61 and Theorems 
15.151 - 15. 16|) . However, until now we considered t as fixed and focused our 
attention on the multiplicative representation (|5.12l - (|5.15|) of the solution 
considered as a function of x. Now we emphasize the dependence of the 
local solution on t, while concentrating on the non-resonant case. Given 
t G C" (such that the coefficients Q p (t),p = 1, . . . , n, are defined for this t), 
and assuming that the matrix Qj(t) is non-resonant for certain j, the local 
solution Yj(x,t) of the equation (|6.1|) corresponding to the singular point tj 
is sought in the form 

(6.3) Yj(x,t) = Hj(x,t)(x - tj) Q *®, x G D tjtP \ {tj}. 

The differentiation rule (|5,8|) for (x — tj)** leads to the differential equation 
for Hj(x, t): 

dHj ( X ,t)) = fQ,{t)H^x,t)-H j{ x,t)Q j{ t) \ 
dx \ x — tj j j j 

where 

Q P (t) 



(6.5) %M)= E 



l<p<n, 



x — t 



r 



the function $>j(x,t) is holomorphic with respect to x at x = tj. The equa- 
tion (|6.4|) does not depend on the choice of branch of the (multivalued, in 
general) function (x — tj)® j ^ . The right hand side of this equation, consid- 
ered as a function of x, has the singularity at the point x = tj. Nevertheless, 
it can be shown that, if the matrix Qj(t) is non-resonant, then, under the 
normalization condition 

(6.6) Hj(x,t)i x=t = I, 
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the holomorphic solution Hj(x,t) of the equation (|6.4j) exists for x from a 
disk centered at tj, and such a solution is unique. More precise formulation 
of this result will be done below. 

In fact, here we repeat what we have already formulated in the Proposition 
15.51 and expressed in the formula (|5.12j) . The only difference between the 
formulae 1)6.3(1 and (|5.12|) is that the notation in ()6.3() reflects the dependence 
on t explicitly. 

Definition 6.2. Let t = (ti,. .. ,t n ) G C™ be given, and assume that the 
matrix Qj{t) is non-resonant. 

The k x k matrix function Yj(x,t), possessing the properties 

i. Yj(x,t) satisfies the differential equation (|6.1|) with respect to x in a 
punctured neighborhood D tj . p \ {tj} of the point tj, (D tjjP = {z£ 
C : \z — tj\ < p} , p is a positive number); 

ii. Yj(x,t) admits the factorization of the form ()6.3J) . where the fac- 
tor Hj(x,t) is a matrix function holomorphic with respect to x in 
the whole (non-punctured) disk Dt j)P and satisfying the normalizing 
condition (|6.6|) . 

is said to be the normalized local solution of the differential equation (|6.1|) . 
corresponding to the singular point tj. 

Proposition 6.3. Let Q p (t),p = 1, . . . , n, be k x k matrix functions, which 
are holomorphic with respect to t = (t\, . . . ,t n ) for t from an open set 
T>, T> G C™. Let t° = (i?, . . . ,t°) be a point from T>, and let 

(6.7) p q (t°) = l min \t° q - t% 1 < q < n. 

Assume that, for certain j, the matrix Qj(t°) is non-resonant. Then there 
exists 5,5 > 0, such that for every t = (t\, . . . ,t n ) from the polydisk 10 
D(t°,<5), the normalized local solution Yj(x,t) of the equation \6. 1]) . which 
corresponds to the singular point tj, exists for x from the punctured disk 
11 -D tj)Pj (tO) \ {tj}, and the factor Hj(x,t) in Kb. 5|) is a matrix function 
holomorphic with respect to x,t for x G D t . p .^ t o^, t G D(£°,<5). Moreover, 
the matrix Hj(x,t) is invertible for these x and t. 

The proof of this proposition will be given in Appendix. 

Remark 6.4. Proposition 16.31 can be considered as a version of Proposition 
15.51 with the emphasis being placed on dependence on t. 

Remark 6.5. Proposition 16.31 belongs to the class of statements of the fol- 
lowing nature: if coefficients of a differential equations and a normalizing 
condition depend on some parameters analytically, then the solution of this 
equation depends on these parameters analytically as well. However, in 

10 The polydisk D(t°,<5) was defined in f^J . 
11 -D tj ,P J ( t 0) = {^eC:| Z -t J |<p J (t )}. 
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Proposition 16.31 the normalizing condition (which can be presented in the 
form ]im x -+tjYj(x,t)(x — tj) - ®^*' = I) is posed at the singular point tj. 
Therefore, Proposition 16.31 is not a consequence of standard general results 
about analytic dependence of a solution of a differential equation on pa- 
rameters. The proof makes explicit use of the non-resonance of the matrix 
Qj(t°) and without this assumption the result fails. This shows that the 
result is rather delicate. 

In section the notion of connection coefficients was introduced. The 
connection coefficients relate the solution of the Fuchsian equation, which 
is normalized at the distinguished point xq, to its local solutions, which are 
normalized at the singular points of the equation. However, in section \5\ we 
did not care about the dependence of connection coefficients on t. Now we 
focus our attention on this dependence, concentrating on the non-resonant 
case. 

Given 12 t°,t° £ T>, and given j, 1 < j < n, let us choose Pj(t°) according 
to (|6.7j) . Then choose 5 > 0, satisfying the condition 

(6.8) 25 < Pj (t°), 

and choose and fix the point Xj G C, satisfying the condition 

(6.9) 8< \xj-t Q j\ <Pj(t°) -S. 




j 




ipj(*°) 








Figure 5. a: disk D t o p .^ t oy 



Figure 5. b: disk D 



t j>Pj (to)- 



With such a choice 
(6.10) Xj G D tjjP]{t0) for every tj G D t o s . 

(See Figure|SJa, where the disk D.o $ is plotted in the light gray color, and the 
dashed circle is the boundary of the disk D t o p .^t ))- Choose a simple path 

1 9 

1? is the open set in where the coefficients Q p , 1 < p < n, are denned and holomorphic. 
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ctj, leading from the distinguished point xq to the point Xj and separated 
from the points t®,p = 1, . . . , t n : 

(6.11) \z — tp\ > 5, Vz£aj, p = l,...,n. 

Then fix the choice of the path ctj. (In particular, let the path ctj be 
independent oft.) From (|6.11j) it follows that for 13 t = (ti,...,t n ) G 
D(t°,S), the condition \x — t p \ > Vx G ctj,p = l,...,n, holds. Thus, 

the function - is holomorphic with respect to x,t for x from a 

l<p<n °° p 

neighborhood of the path ctj and for t G D(t°,5). Therefore, the solution 
Y(x,t) of the differential equation (|6.1j) that satisfies the initial condition 

(6.12) Y(x,t) x=X0 =I 

can be continued analytically along the path ctj from a neighborhood V xo 
of the initial point xq to a neighborhood V Xj of the point Xj. (V Xj does 
not depend on t £ D(t , <5)). We denote by Y aj (x,t) the function in V Xj 
which is defined by means of such an analytic continuation. Since the initial 
condition in the (regular) point xq and the path ctj do not depend on t, then 
for every x G ctj, the solution Y a .(x,t) is holomorphic with respect to t for 
x G V Xjj , t G D(t°,S). In particular, the matrix Y a . (xj, t), as well as the 
inverse matrix Y~^{xj,t), is holomorphic with respect to t for t G D(t°,8). 
It is clear that the solutions Y a i (x,t) and Y a >>{x,t) corresponding to the 
analytic continuations into V Xj along two different paths ctj and a'j (with 
the same endpoints xq and a%), are related by the equality 

(6.13) Y a ,,(x,t) = Y a ,(x,t)-M 7 (t), 

where M 7 (t) is the monodromy matrix for the equation ()6.1|) corresponding 

to the loop 7, 7 = (a^) a!-, composed from the paths ctj, a". (The loop 7 
does not depend on t, but the equation depends. So, the monodromy matrix 
Mj(t) can depend on t). 

Now we turn our attention to the local solution Yj{x,t), corresponding 
to the singular point tj. Assume that the matrix Qj(t ) is non-resonant. 
Choosing a smaller 5, if necessary, we can apply Proposition 16.31 According 
to this proposition, the local solution Yj(x, t) exists in the punctured disk 
D t p (joy \ {tj}. In particular, it is defined at the point Xj, and even in the 
neighborhood V Xj , if this neighborhood is small enough. The local solution 
Yj(x,t) is the product of two factors: Hj(x,t) and (x — tj)®^*'. Under 
the normalizing condition (|6,6[) . the first factor is determined uniquely. The 
second factor is not quite unique. In general, the function (x — tj)® j ^ is a 
multivalued function of x. This function, and hence the solution Yj(x, t), are 
determined up to right factor exp{2iriQj(t)}. To avoid this non-uniqueness, 
we must come to an agreement on the value of arg(x — tj) for some fixed 



13 D(i°,<5) = D f0 g X • • ■ X D t o j is a polydisk. 
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x, for example, for x = Xj. However, now t = (ti, . . . , t n ) is not a constant 
but a variable with values in the polydisk D(t°,5). In particular, the j-th 
coordinate tj takes the values in the disk D t o p . In order for the function 

(xj —tj)^ j ^ to be holomorphic with respect to t, we need to choose arg(xj — 
tj) coherently with respect to tj. 

Let us choose somehow and fix the value arg(xj — tj), say 

(6.14) arg(ay -t°j) = $°j. 

For tj G D t o s , choose the value of arg(xj — tj) according to the rule: 

(6.15) i? - < axg(xj - tj) <$ + vr/2. 

The choice (|6,15[) is possible and unique (see FigureEla). With such a choice, 
the matrix (xj — tj) 

Q 3 (t) ig 

determined uniquely and is holomorphic with 
respect to t for t £ D(t°, 5). Therefore, the matrix Yj(xj,t) = Hj(x, t)(xj — 
^.^Qj{t) - g detefjixffled uniquely and is holomorphic with respect to t for 
teD(t°,6). 

Moreover, the function Yj(x,t) — the local solution corresponding to the 
singular point tj — is determined uniquely on the universal covering 
cov(D tj p .(to^ tj j, Xj) with the distinguished point Xj. The conditions ()6.14j) - 
()6.15|) . together with the condition ()6.6|) . can be considered as normalizing 
condition for the local solution Yj(x, t): the condition (|6.6() normalized the 
factor Hj(x,t)), the condition ()6.15|) normalized the factor (x — tj)®'^ 
coherently with respect to tj. 

Let us define the connection coefficient C aj (t) as 

(6.16) C a . (t) = Y-\xj,t) ■ Y Xh<Xj ( Xj ,t) 

Since both factors in the right hand side of (|6.16j) are holomorphic with 
respect to t £ D(t°,<5) and invertible, the matrix C aj (t) is holomorphic 
with respect to t E D(£°,<5) and invertible as well. According to ()6.13|) . the 
connection coefficients C' a .(t) and C'^.(t), corresponding to two different 
paths a'j and a", are related by the equality 

(6.17) C a »(x,t) = C a >.(x,t) ■ M 7 (t), 

where M 7 (t) is the monodromy matrix for the equation (|6.1|) corresponding 
to the loop 7, 7 = (a'j) a". 

Since the matrix C aj (t) does not depend on x, the product Yj(x,t) ■ C a . 
is a solution of the differential equation (|6,1|) for x G D t in particular, 

for x E V Xj . The function Y x . tCt . (x, t) is also a solution of the differential 
equation (|6.1|) for x € V Xj . Both these solutions coincide at the point Xj. 
Hence they coincide for all x G V X] . Therefore, the product Yj(x,t) ■ C aj 
can be considered as the analytic continuation of the solution Y x . fa . (x,t) 
from V X] to the punctured disk D t . p .^ \ {tj} (or, better to say, to the 
universal covering cov(D t .^ p .^ t o- ) \ {tj};xj) of this punctured disk with the 



RATIONAL SOLUTIONS OF THE SCHLESINGER SYSTEM 



33 



distinguished point Xj). In its turn, the solution Y Xjt0l .(x,t), which is de- 
fined in V x , was obtained by the described above analytic continuation of 
the solution, "living" in V Xo and satisfying the initial condition (|6.12j) . In 
other words, the solution Yj(x, t) ■ C a . can be considered as a solution which 
is obtained as the result of the "two step" analytic continuation from V xo 
to cov(D t . pj ( t o} \ {tj};xj) of the solution, satisfying the initial condition 
(|6.12j) . However, it is better to talk not about the two-step analytic contin- 
uation but about the analytic continuation along the composed path. The 
representative of a point of the universal covering cov(D t ^ p .^ t o- ) \ {tj}; Xj) is 
a pair where x is a point of the punctured disk D t ^ p .^ \ {tj}, and 

fi C {D t . p .^ \ {tj}) is a path starting at the point Xj and ending at x. The 
path corresponding to the above mentioned two-step analytic continuation 
is the composition 14 a = fi ■ aj. Thus, to describe the dependence on t 
of the solution Y aj (x,t), which is normalized at the initial point xq by the 
condition 

(6.18) Y aj (x,t) ]x=XQ =I, 

in the neighborhood of the singular point tj, we continue analytically this 
solution from the neighborhood V xo of the initial point to the punctured 
neighborhood D t .^ p .^ t o^ along those paths a which reach first the point Xj 
(xj does not depend on tj and is common for all these paths), and then 
go, remaining within the punctured neighborhood D t ^ p ^ t o^ \ {tj}, to the 
destination points . 

Definition 6.6. Such a family 15 {Y aj . (x, *)}t g D(t ,<5) °f analytic continua- 
tions (with respect to x) of the normalized by (|6.18j) solutions of the equa- 
tion (J6.1|) from a neighborhood V Xo of the initial point xq to the fam- 
ily {D t . p ^ t o^}t J ^D t0 s of the punctured disk neighborhoods of the points 

tj,tj € D t o s , is said to be coherent with respect to t. The path ay along 
which such continuations are carried out is said to be the path determining 
the coherent family of analytic continuations. 

The solution 

(6.19) Yj(x,t) ■ C aj {t) = H 3 (x,t)(x - tj) Q ^C a] {t) 
can be presented in another form - more suitable for our goal: 

(6.20) Y aj (x,t) = H aj (x,t)(x - tjf^W, 

where 
(6.21) 

a). A aj (t) = C-j{t)Qj{t)C aj {t); b). H aj {x,t) = H 3 {x,t)C a] (t). 



The path a starts at the point xq and ends at the point x. Traversing the path a, we traverse 
first the path ctj, and then the path fi. See Figure l5lb 
^Indexed by the parameter t £ D(t°,S) 
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Of course, the relation (|6.17JI leads immediately to the relations 
(6.22) 

a). A a ,,(t) = (M y (t))- 1 A a/ (t)M y (t), b). H a ,,(x,t) = H^x^M^t), 

where a'j and a" are two paths with the same endpoints xo and Xj, and M 7 
is the monodromy matrix for the equation (|6,1|) corresponding to the loop 

7, 7 = («j)~ lQ i- 

Definition 6.7. The matrix functions H aj (x,t) and {x — tj) Aa ^ t " > are said 
to be the regular factor and the principal factor of the representation (|6.20f) 
of the solution Y a .{x,t) near the point tj. 

Let us summarize the results, obtained above, as 

Theorem 6.8. Let the k x k matrix functions Q p (t),p = l,...,n, which 
appear in the Fuchsian system (|6.1I) be holomorphic with respect to t = 
(ti, . . . , t n ) for t from an open set T>,T> G C". Let t° = (t®, . . . , t^) be a point 
from T>. Assume that, for certain j, the matrix Qj(t°) is non-resonant, and 
let Pj(t°) be defined according to \6.T\ l) (with q = j). 

Then there exists a number 5, < 5 < ^pj(t°), such that, fort € D(t°,S), 
the solutions {Y aj (x, t)}t£D(t°,6) of the holomorphic family of equation \6.1\) 
which are distinguished by the properties: 

i. Each of these solutions is normalized by (|6.18|) at the initial point xq 

ii Each of these solution is continued analytically from a neighborhood 
of xq to the family {D t . ( t o)}t,eD - of the punctured disk neighbor- 

hoods of the points tj,tj G D t q and such analytic continuations 

are performed coherently with respect to t G t G D(t°,5) , 

are representable in D t . p ^ t o^ \ {tj} in the form (|6,2U|) . where: 

(1) The matrix function H aj (x,t) — the regular factor of the represen- 
tation (|6.2U|) — is holomorphic together with its inverse H~j(x,t) 
with respect to x,t for x G D t ^ p ^ t o^,, t G D(t°,S); 

(2) The exponent A a . (t) of the principal factor (x — tj) Aa j ^ is a matrix 
function holomorphic with respect to t from the polydisk t G D(i°,<5); 

(3) The exponent A a .(t) is similar to the matrix Qj(t), §b\21i a). and 
the similarity matrix C a .(t) can be chosen 16 to be holomorphic for 
t G D(t°,<5) . 

(4) By a different choice of the paths a'j and a", along which the coher- 
ent analytic continuations are carried out, the matrices A aj {t) and 
H a \x,t) are transformed according to the rule (|6.22|) . where M^{t) 



The similarity matrix is not unique: it can be multiplied by an arbitrary matrix commuting 
with Qj(t). 
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is the monodromy matrix corresponding to a loop 7, constructed 
from the paths a'j and a". 

Let ctj be a path, separated (as in (|fj.lljl ) from the singular points t p , 1 < 
p < n, and leading from the point xq to a point Xj E (C\ ^tj.p^t )) \ ^t°-,Si 
where the intersection is taken over all tj G D t o$. This path can be used 

for the coherent analytic continuations of the normalized by ()6.18j) solution 
to the neighborhoods of the point tj G D t o s . The path cxj also generates 
a loop 7 Q . which goes around the disk neighborhood D t o s. This loop does 

not depend on t and, for every t = (ti, . . . ,t n ), t G D(t°,5), represents an 
element of the fundamental group 7r(C \{t±, . . . , t n }; xq). (Compare to the 
Definition I5,12|) , The loop j aj is the composition of the path ctj, of the 
circular arc \i, starting and ending at the point xj, and of the path aj 1 (the 
path aj 1 is the path aj, traversed in the "opposite direction"). In more 
detail, /j, is the circular arc of the radius \xj — t®\, centered at the point tp 
which starts at the point Xj, makes the full counterclockwise turn around 
the point tp and ends at the same point Xj (or, better to say, at the point 
i° + e 2 ™{ Xj - f°)). Since Xj G {(] D tjtPj (tP) ) \ D t o >5 , n C D thp . [lP) \ {tj} for 
every tj G D t o s . Since in \ {tj} the representation (|6.20j) holds, 

Theorem I5.18lf can be formulated in the following "t-dependent" version: 

Theorem 6.9. Assume that the matrix Qj(t°) is non-resonant. Let Oj 
be a path leading from the initial point xq to a punctured neighborhood of 
the singular point tp and let j aj be the loop generated by the path aj (as 
described above). 

If t G D(t°,5), and 5 > is small enough, then the monodromy matrix 
M 7a .(t) for the equation corresponding to the loop j aj , and the expo- 

nent A a . (£) of the principal factor of the representation tyS.'l{)\ of its solution 
Y aj (x,t) are related by the equality 

(6.23) M la] (t) = exp{2mA aj (t)}. 

The question arises, whether the monodromy matrix M 7q . (t) determines 
the exponent A aj (t) uniquely? In particular, does the equality A a p{t') = 
A a .(t") follow from the equality M ltx .{t') = M^(t"), ift',t" G D(t°,S), 
and 5 G is small enough? 

The latter question is closely related to the description of the set of solu- 
tions of the matrix equation (with respect to ^4): 

(6.24) exp{2iriA} = M, 

where M is a given k x k matrix. A description of the set of solutions of 
(|fi.24j) can be found in |(irantj . Chapter VIII, §8. The invertibility of M is 

^ Both paths, oi- and a", have the common starting point xq. If their endpoints coincide, then 
the loop 7 is just the composition 7 = (a^) — 1 a". In the general case, in which these endpoints 
may be different, 7 = where v is a path, connecting the endpoints. 
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the necessary and sufficient condition for the solvability of (|6.24j) . Given 
M, det M 7^ 0, the set of solutions of 1)6. 24|) is always infinite. Our prime 
interest is in conditions under which a certain solution Aq of the equation 
(|6.24j) is isolated. In general, non-isolated solutions may exist. For example, 
let M = diag[l,l] , Aq = diag[0, 1], and let T be an arbitrary invertible 
matrix. Then exp{27ri(T~ 1 AqT)} = M. If T is non-diagonal matrix, then 
T~ 1 AqT Aq, but if T is close to /, then T~ 1 AqT is close to Aq. 

Lemma 6.10. Given a non-resonant matrix Aq, there exists e > such that 
from the conditions exp{2iriAQ} = exp{2iriA} and \\A — Aq\\ < e it follows 
that A = A . 

PROOF. Since Aq is non-resonant, there exists an open (not necessarily con- 
nected) set O, O C C containing the spectrum of Aq such that the map- 
ping £ i— > exp{27ri£} is univalent (schlicht) in O. If e is small enough, the 
set O also contains the spectra of all matrices A satisfying the conditions 
ll-A— Aj|| < e. The set Q = exp{2iriO} is an open set containing all the spec- 
tra of the matrices exp{27riA} with 1 1 ^4 — ^4o 1 1 < £ i hi particular, the spectrum 

of the matrix exp{27riylo}- If Lnw, Ln : Q *— > O, is a function, inverse 

2iri 2ixi 

to the function exp{2iri(^} : O i— > Q, that is Ln(exp{27r«C}) = Q V£ E O, 

^ 2iri 

then Ln(exp{2-7ri^4}) = A for all matrices A whose spectra are contained 

2iri 

in O, in particular, V A : \\A — Aq\\ < e. Q.E.D. 



Since the matrix A aj (t°) is similar to the matrix Qj(t°), the following result 
is the direct consequence of Lemma 16.101 and Theorem 16.91 

Theorem 6.11. Let the matrix Qj(t°) be non-resonant, and let A aj (t) and 
Mj a . (t) be the same as in the formulation of Theorem \b'.!A If 5 > is small 

enough, and if for some t',t" € D(t°,S), the equalities M loi .{t r ) = M la , (t") 
hold, then the equalities A a .(t') = A a .(t") hold as well. In particular, if 
M la .(t) = M Ja .(t°) in D (*°^); then A a .(t) = A aj (t°) in D(i°,<5) as well. 

Remark 6.12. If the matrix Qj(t°) is resonant, then, trying to formulate a 
"t-dependent" version of Theorem 15.161 we should seek the representation 
of the solution Y a . (£) near £° in the form 
(6.25) 

Y a .{x,t) = H aj (x,t)(x - tj) Za i {t \x - tj) Aa i {t \ x G cov (D tj:P \ {tj}, Xj ), 

where the matrices H a . (x, t), H~j(x, t) are holomorphic for x G D t p ^ 
and t G D(i°, pj(t )), the matrices Z a . (£), A a . (t) are holomorphic for t G 
D(£°, pj(t )), and moreover the matrices Z a .(t) are diagonalizable , the 
eigenvalues of Z a .{€) are integer, the eigenvalues of A aj (t) and the eigen- 
values of Qj(t) are related by equalities analogous to the equalities (|5.29|) . 
However, if the matrix Qj(t°) is resonant, the representation of the form 
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(|o.25j) . with holomorphic Z a . (£), A aj (i), does not exist in general. The rea- 
son is that it is possible for the matrix Qj(t°) to be resonant, while all the 
matrices Qj(t) for t ^ t are non-resonant. (The property of resonance is 
unstable under small perturbations). According to Theorem 15.151 the solu- 
tion Y aj (t) must be of the form (j(12()j) for t ^ t . According to Theorem 
15.161 it must be of the form (|6.25j) for t = t°, where in general Z a .(tP) ^ 0. 
However, such a bifurcation is incompatible with the holomorphy of Z aj (t). 

7. ISOPRINCIPAL AND ISOMONODROMIC DEFORMATIONS. 

Consider a Fuchsian equation (|6.1|) . where k x k matrices Q p (t), 1 < p < 
n, are holomorphic in t for t from a domain T>, T> C C™. The equation 
H6.1J) is considered as a differential equation with respect to the variable x, 
and t = (ti, . . . ,t n ) is considered as a parameter. Recall that we use the 
terminology "the deformation of Fuchsian equation" for a Fuchsian equation 
depending on a parameter 18 . 

Definition 7.1. Let (|6.1[) be a deformation of Fuchsian equation, where 
k x k matrices Q p (t), 1 < p < n, are holomorphic in a domain T>,T> C C™, 
and let t° = (i§, . . . be a point of X>. The deformation is said to be 
isomonodromic with respect to a distinguished point xq at the pole loci t° 
if for every loop 7 with the distinguished point xq, 7 C C \ {t^, . . . , t^}, 
there exists S > 0, 5 = 5(t°,j), such that the monodromy matrix M~(t) 
of the deformation (jfi.lj) along the loop 7 does not depend on t in D(t , 5): 
Mj(t) = M 7 (t°) Vt G D(t°, 5). (Of course, we assume that 5 < dist(t°,7), so 
that 7 C C \ {ti, . . . , t n } for t G D(t°, <5), and the monodromy matrix M 7 (t) 
is well defined for such t). The deformation is said to be isomonodromic with 
respect to xo in T> if it is isomonodromic with respect to xq at every pole 
loci t eT>. 

Remark 7.2. To check that the deformation (|fi.l|) is isomonodromic at i°, 
there is no need to examine the monodromy matrix M 7 (t) for all loops 7. It 
is enough to choose loops 7^, j = 1, . . . , n, 7^ C C \ {t^, . . . , t° n }, generating 
19 the fundamental group 7r(C \ {t^, . . . , t^}, xq), and to check the the mon- 
odromy matrices M 7j . (i) do not depend on t for £ G D(£°, <5), <5 > is small 
enough. 

Remark 7.3. If 7 is a loop with the distinguished point xo, and tj 7,J = 
l,...,n, then the monodromy matrix M 7 (.) is holomorphic at the point 
t = (ti, . . . , t n ). From the uniqueness theorem for holomorphic functions and 
from elementary topological considerations it follows that if the deformation 
(|6.1j) is isomonodromic at some t G T>, then it is isomonodromic at every 
t G T>. 

18 We consider only such deformations, where the loci t = (ti, ■ ■ ■ , t n ) serve as a parameter, 
and the dependence on the parameter is holomorphic. 

l^If t £ D(t°,5) and <5 > is small enough, then such loops fj satisfy the condition 7,- C 
C \ {t± , . . . , tn} and generate the fundamental group 7r(C \ {t\ , . . . , t n }, xq). See Figure FT1 
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Definition 7.4. Deformation (jfi.lj) is said to be isospectral, if for every 
j, j = 1, ... ,n, the spectrum of the matrix Q./(t) does not depend on t in 
T>: spec Qj(t') = spec Qj(t") (the sets coincides "with multiplicities") for 
every t',t" G X>. 

Theorem 7.5. // a deformation (|6.1j) is isomonodromic with respect to a 
distinguished point, then it is also isospectral. 

PROOF. For given j, the coefficients of the characteristic polynomials of 
the matrix Qj(t) are holomorphic with respect to t in X>. Therefore, it is 
enough to prove that for some 5 > 0, this polynomial does not depend on 
t for t G D(i°,<5) or, what is the same, that the spectrum specQj(t) does 
not depend on t for t S D(i°,<5). Let jj be a loop with the distinguished 
point xo wich makes one turn counterclockwise around the singular point 
t® and makes no turns around other singular points t®, p = l,...,n,p^j, 
as it was described in the formulation of Theorem 15.191 The monodromy 
matrix M 7j (i) does not depend on t for t which are close to t° (In general, 
M 7j .(t) = M 7i (t°) if it is possible to pass from t° to t in "D so that no t p 
crosses the loop 7j). According to Theorem 15.191 the set exp{ spec Qj(t)} 
does not depend on t in D(i°,<5). However, on the one hand, the spectrum 
specQj(i) depends on t continuously, and on the other hand, the mapping 
C i— > exp{C} is locally univalent (schlicht). Therefore, the set spec Qj(i) does 
not depend on t in D(t°,5) if 5 is small enough. Q.E.D. 

Let (|6.1|) be an isomodromic with respect to deformation of Fuchsian 
equation, and let for a certain j, the matrix Qj(t°) be non-resonant 20 . 
Let ay, ay C C \ {tj 1 ; • • • be a path leading from the point xq to a 
small neighborhood of the point t®. According to Theorem I6.11( for t 6 
D(i°,<5), the solution Y aj (x,t), normalized at the distinguished point xo by 
the condition (|6,12|) and continued analytically to the neighborhood of t®, 
is representable near L- in the form 

(7.1) Y a .(x,t) = H aj (x,t)(x-t J ) Aa i, xeD tj , p ,teD(t°,5), 

where 5 and p are some positive numbers, the regular factor H(x, t) is holo- 
morphic in x,t and invertible for x £ ^tj.p-, t & D(t°,<5), and the exponent 
A aj of the principal factor (x — tj) a i does not depend on t. Thus it is 
natural to give the following definition: 

Definition 7.6. The deformation (|6.1j) is said to be isoprincipal with respect 
to a distinguished point xq in the narrow sense at the pole loci t° = (t®, . . . , £^), 
(t° € T>), if for every j = 1, . . . , n and for some paths a±, . . . ,a n leading 
from xo to neighborhoods of the points t\, . . . ,t° respectively, the solution 
Y a .(x,t) normalized at xq by (|6.12|) is representable in the form (|7.1|) near 



^0 According to Theorem 17,51 the matrix Qj (t) in non-resonant for every t£P. 
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tj, where the matrices A aj do not depend on t for t £ D(t , 5), 5 > is 
small enough. 

We emphasize that we do not require in this definition that the matrices 
Qj(t°) are non-resonant. Of course, in order for this definition to be correct, 
we must check that it does not depend on the choice of the paths ay. This 
can be done easily. Below (Lemma I7.13[) we prove this independence in a 
more general situation. Let us summarize the results in the form of two 
lemmas: 

Lemma 7.7. In the non-resonant case (i.e., if all the matrices Qj(t°) are 
non-resonant), every isomonodromic deformation is isoprincipal in the nar- 
row sense as well. 

This lemma is the immediate consequence of Theorem 16.111 The converse 
is true always, without any assumption of non-resonance. 

Lemma 7.8. If the deformation (|6.1|) is isoprincipal (in the "narrow sense") 
at some t° (no assumptions on the matrices Qj(t°),j = 1, . . . ,n are made), 
then this deformation is isomonodromic as well. 

PROOF. If the solution Y a . (x,t) of (|6.1j) is representable in the form (|7.1|) 
with a constant A aj , then M 7a .(t) = exp {27riA aj }, where M 7a . is the 
monodromy matrix along the loop "y a . generated by the path ctj. Therefore, 
M 7 does not depend on t. Q.E.D. 

Remark 7.9. If a deformation (|6.1|) is isomonodromic, but some of the ma- 
trices Qj{t) are not non-resonant, then this deformation may be not isoprin- 
cipal. As we shall see later, there is a rich class of deformations of Fuchsian 
equations all fundamental solutions Y(x, t) of which are rational functions 
of x. Thus, such deformations are isomonodromic. (Their monodromies are 
trivial). However, most of these deformations are not isoprincipal. 

From Lemmas 17.71 and !7.8l it follows that in the non-resonant case the classes 
of isomonodromic deformations and of isoprincipal (in the "narrow sense" ) 
deformations coincide. This is a useful observation. On the one hand, it 
is the property of " isoprincipalness" that implies the Schlesinger system for 
the matrix functions Q p (t) appearing in the deformation. (We prove this 
implication below). On the other hand, there are methods, which are based 
on solving the so called Riemann-Hilbert problem, that allow to construct 
deformations of Fuchsian equations with the prescribed monodromy, but 
not with the prescribed principal factors of the solutions Y a . . In particular, 
such methods allow (under non-resonance condition and certain other re- 
strictions) to construct directly isomonodromic (but not isoprincipal) defor- 
mations. The proof of the fact that the resulting Qj(t) satisfy the Schlesinger 
system uses essentially the equivalency of these two classes of deformations 
in the non-resonant case. 

However, the above given "narrow" definition of the isoprincipalness, 
which is quite natural in the non-resonant case, seems to be a little bit 
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artificial and too restrictive in the general case. Indeed, in the general case 
the solutions Y aj (t) of the equation (|6.1j) are of the form 1)6. 25|) . where expo- 
nents Z aj (i) and A aj (t) depend on t in general and posses some additional 
properties. In particular, the matrix A a Jt) is non-resonant, the matrix 
Z aj (£) is diagonalizable, with integer eigenvalues. (These exponents may be 
non-holomorphic, and even discontinuous functions of t. See Remark 16. 12|) . 
If these exponents are constant then we have good reason to call such a 
deformation isoprincipal. However, we have no need to go into detail. All 
that is important for us is that the principal factor (x — tj) Za j ■ (x — tj) Aa i 
depends on the difference x — tj only if the exponents Z aj and A aj do not 
depend on t. 

Definition 7.10. The deformation 1)6.1)) is said to be isoprincipal with respect 
to a distinguished point xq (in the wide sense) at a pole loci t° = (t®, . . . , 
(t° € T>), if for every j = 1, . . . , n and for some paths a\, . . . , a n leading 
from xo to neighborhoods of the points t-jj, ... ,t° respectively, the solution 
Y aj {x,t) normalized at xo by 1)6.12)1 is representable in the form 

(7.2) Y aj (x,t) =H a .(x,t)E a .(x-t j ), x E D tjtP \ {tj}, t £ D(t°, 8), 

where the matrix functions H aj (x,t), (H aj (x,t)} 1 are holomorphic in 
x,t for x E Dt JjP ,t € D(t°,5), (p and 5 are some positive numbers), and 
E aj {0-,3 = 1, ... ,n, are matrix functions possessing the properties: 

(1) E aj (Q are holomorphic on the Riemann surface of logarithm 21 and 
take invertible values there; 

(2) E aj (() do not depend on t and x. 

Remark 7.11. Of course, to formulate Definition I7.1UI more accurately, we 
have to describe the geometric configuration more accurately. For every 
j, we have to choose the point Xj according to (j6.9|) . assuming that 5 is 
small enough, in particular, (|6.8j) is satisfied. (The value Pj(t°) is defined 
according to (|6.7|) ). By such a choice of Xj, the condition (|6.10[) is satisfied. 
Then we have to choose the arguments &ig(xj — t®) somehow, (jti.14)) . and 
then &icg(xj — tj) for tj € D t o s according to ()6.15j) . A point of the universal 
covering cav{D t . p .^ t a^ \ {tj}, Xj) is a pair (x, p), where x is a point from the 
punctured disk D t . p .^ t o^) \ {tj}, and p is a path which starts at the point Xj 
and ends at the point x, p £ D tj)Pj (j.o^) \ {tj}. The path p, together with the 
choice of &ig(xj— tj), determines the value of arg(x— tj), and hence, the value 
of E aj (x — tj). On the other hand, the pair (x, p) G cov(D t .^^ t o- ) \ {tj}, Xj) 
determines the value of the analytic continuation Y a . along the path ctj at 
this point (x,p). By definition, the value Y a .{(x, p),t) is the result of the 
analytic continuation of the solution normalized at x$ by (|6.18|) along the 
composed path p ■ aj from a neighborhood of xq to the point x. Since both 



The Riemann surface of logarithm is the universal covering cov(C \ {0}) . Considered as 
functions defined in C \ {0}, E aj (£ ) are multivalued in general. 
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values, Y a . (x, t) and E aj (x — tj) are determined for every choice of the path 
fi (these values may depend on fx), the value H a .{x,t) = Y a . (x, t) [E a . (x — 
tj)) 1 is determined for x £ F> t ^ p .^ t o^ \ {tj}, and it may depend on fx. The 
definition I7.1UI means firstly, that the value H a . [x, t) is determined only 
by the point x S D t .^ p .^p^ \ {tj} itself, but does not depend on the path 
fx leading to this point, thus the function H Qj (x, t) as a function of x is 
single-valued in the punctured disk D tjtPj ^ \ {tj}. Secondly, this definition 
means that the singularity x = tj of the function H aj {x,t) is removable, 

that is there exists the limit lim x ^^ )X ^f ,H aj {x,t) = H aj (tj,t). Thirdly, 
the definition I7.1UI means that the value H a .{tj,t) is an invertible matrix. 
(For x ^ tj, the values H aj (tj,t) are invertible, because both the values 
Y aj (x,t) and E a .(x — tj) are invertible). So, the accurate definition of the 
isoprincipal deformation is rather bulky. 

Below ( Lemma 17.13(1 we prove that the notion of isoprincipal deformation 
is well defined, that is the property of the deformation to be isoprincipal at 
some point does not depend on the choice of the paths ctj. 

Theorem 7.12. If (jfi.lj) is a deformation of Fuchsian equation which is 
isoprincipal with respect to a distinguished point in the wide sense (i.e, in 
the sense of Definition \ 7. at a pole loci t° = (t®, . . . ,t^), (t° $LT>), then 
the matrices E aj ((^),j = l,...,n, representing the principle factors, must 
satisfy the monodromic relations of the form 

(7.3) E a . (Ce 2 ™) = E aj (C)M Qj , V( + 0, 

where M a . are constant (not depending on Q) invertible matrices. For the 
monodromy matrix M 7q . (£) of the equation (|6.1j) along the loop r y aj gener- 
ated by the paths ctj, (see Figure^), the equality 

(7.4) M la ,{t) = M a ., t€D(f,S), j = l,...,n. 

holds. In particular, this deformation is isomonodromic with respect to xq. 

PROOF. For fixed t € D(i°,5), and for x € Dt tP \ {tj}, consider two ma- 
trix functions of x: Y aj (x,t) and Y aj {tj + (x — tj)e 2m ,t). Here x and 
tj + {x — tj)e 2m are considered as points running over the universal cov- 
ering cov(D tj! p \ {tj } , Xj ) of the punctured disk D tjtP \ {tj} with the dis- 
tinguished point Xj. These two points of the universal covering have the 
same projections in the complex plane. Therefore, the matrices Y aj {x,t) 
and Y a .(tj + {x- tj)e 27ri ,t) can be considered as the values of two differ- 
ent solutions of the same linear differential equation at the same point. 
From &etY{x,t)\ x=x . ^ (see (02]) ) it follows that det Y a . (x, t) ^ 0, 
detY aj (tj + (x — tj)e 2m ,t) ^ 0. Hence, these solutions are "proportional" 
as functions of x: 
(7.5) 

Y aj {tj + {x- t 3 )e 2m ,t) = Y aj {x,t)M aj {t), x G D t . jP \ {tj},t G D{t°,p), 
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where M a .(t) are constant (with respect to x) non-degenerate matrices, 
which may depend on t. Since the factor H a . is univalued {H aj (tj + {x — 
tj)e 27Ti ,t) = H a .(x,t)), from and (J73J) it follows that 

(7.6) E aj ((x - tj)e 2 n = E a . (x - t 3 )M a . (t) 

Since x is an arbitrary point of D t . p \{tj}, x — tj is an arbitrary point of 
C \ {0}. Since the function E aj does not depend on t, the matrix M aj does 
not depend on t as well. Moreover, the monodromic relation (|7.3|) holds at 
lest for C : < |C| < p. However, the function E a . is holomorphic in the 
whole Riemann surface of logarithm. Therefore, the relation (|7.3j) holds for 
all C : |Cl > 0. 

The expression for the monodromy matrix M 7a . (i) along the loop j aj 
was already obtained (see (|5.19|0 : 

(7.7) M 7aj (t) = V",.,.„ (.r r /.*;• 1 • Y Xita .{tj + (x 3 - tj)e 2 ™,t). 

The equality (|7.4|) is the consequence of (|7.5|) and (|7.7|) . Q.E.D. 

Let for a certain p, a p be a path leading from the point xq to a neighborhood 
of the point t p . Let a' p be another path leading from xq to the same neighbor- 
hood. Let Y ap and Y a i^ be the analytic continuations of the normalized at xq 
solution of (|6.1|) along the paths a p and a' p respectively. Then the solutions 
Y ap and Y a > p are related by the relation Y a r (x,t) = Y ap (x,t) ■ M 7 (t), x € 
Dt ,p \ {tj}, where M 7 (t) is the monodromy matrix along the loop 7, con- 
structed from the paths a p and a' p . (See footnote I17|) . If the deformation 
(|6.1j) is isoprincipal in the sense of the definition 17.101 bv a certain choice of 
the paths aj,j = 1, . . . ,n, then, according to Theorem 17.121 this deforma- 
tion is isomonodromic. In particular, the monodromy matrix M 7 (t) in (|7.8j) 
does not depend on t: M 7 (t) = M 7 , thus 

(7.8) Y a , p (x,t)=Y ap (x,t)-M 7 , xGA P) p\{^}, 
Therefore, the solution Y a > p is representable in the form 

(7.9) Y a , p (x,t) = H a , p (x,t)-E a , p (x-t j ), x <E D tjtP \ {tj},t e D(tj,p), 

with H a / p (x, t) = H ap (x,t)Mj, E a / p (() = M~ 1 E olp (()M~ / . Since the matrices 
EacpiO an d -^7 do not depend on t, the matrix E a ' p (C) does not depend on 
t as well. So, we prove the following 

Lemma 7.13. The notion of the isoprincipal deformation at a point is well 
defined. If the deformation <\6.1^ is isoprincipal at the point t° for some 
choice of the paths aj, j = 1, . . . , n leading from xq to neighborhoods of the 
singularities t®, then this deformation is isoprincipal at this point for any 
other choice of such paths. 
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Definition 7.14. The deformation (|6.1|) . where the coefficients Q p (t),p = 
1, . . . , n, are holomorphic for t £ T>, T> is an open set in C™, is said to be 
isoprincipal (in the wide sense) in X> if it is isoprincipal at every point t E *D. 

Theorem 7.15. If the deformation ()6.1|) is isoprincipal with respect to xq 
at some t°, t° £ T>, and the open set T> is connected, then this deformation 
is isoprincipal with respect to xq in T> (i.e., at every t £T>). 

We omit proof of this theorem. 

8. ISOPRINCIPALNESS OF A DEFORMATION IMPLIES THE 
SCHLESINGER SYSTEM FOR ITS COEFFICIENTS 

Lemma 8.1. Let the deformation 22 ()6.1|) be isoprincipal with respect to a 
distinguished point xq (in the wide sense) at some pole loci t° , t° € T>, 
and let Y(x,t), x £ cov(C \ {t\, . . . ,t n }),t E T>, be the solution of the dif- 
ferential equation H6.1\) which is normalized by the condition (|6.12|) at the 
distinguished point xq. 
Then the equalities 

(8.1) ?Xp±l = Qk{t) {l L-W X| t), k = l,...,n 

dt k \x -t k x-t k J 

are satisfied. In particular, ifxo = oo, then the equations ()4,lb|) are satisfied. 

PROOF. The key to the proof are the (trivial, but nevertheless very useful) 
equalities 

(8.2) A. Ea] {x-t j ) = --^E aj {x-t j ), j = l,...,n, 

which express the fact that the functions E a - (x — tj ) depend only on the 
difference of the arguments x and tj. Expressing the logarithmic derivative 

f rom differential equation 116.11) : Y^ 1 = - 

dx dx x — t 



we 



P "P 



derive its behavior in a neighborhood of the singular point tj\ 
dY 

(8.3) — (x,t)Y-\x,t) = 



Qj(t) 

X — tj 



+ (a function of x, holomorphic with respect to x), 

x e D tj!P , j = l,...n. 



(D tjjP = {z £ C : \z — tj\ < p}, p is a positive number). We restrict 
our consideration to t £ D(t°,<5), where 5 > is small enough, so that 
the representations (|7.2|) of the solution Y{x,t) hold in the appropriate 



22: 



The condition Yl Qp(t) = 0, t £ T>, is assumed. 

Kp<n 
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neighborhoods D tjtP \{tj} of the singular points tj. Using the representation 
(|7.2j) . we see that 

(8.4) 

^L( x ,t)Y- 1 (x,t)=H a] (x,t) (±E aj (x-tj)E-}(x-tj)\ H-i(x,t) 

+ (J-H aj (x,t))\ H-j{x,t). 

Since the function H a .(x,t) is holomorphic and invertible in an entire (non- 
punctured) neighborhood of the singular point tj, comparing (|8,3|) and (|8.4jl . 
we obtain 

(8.5) H a .(x,t) (^E aj {x - tj)(E a .{x - tj))" 1 } H a] (x,t)- 1 = 



Qj 



j 



X tj 



+ (a function of x, holomorphic with respect to x), 

x e D t p , j = l,...n. 



QY 

Let us investigate the logarithmic derivatives — — (x,t)Y ( func- 

dt k 

tions of x. First of all we remark that for fixed t, these function are single 
valued and holomorphic with respect to x in the domain C\ {t±, . . . , t n }. In- 
deed, let O be a simply connected open set in C separated from the points 
t\, ...,t° n : dist(C, {t?, . . .,t° n }) > 0. Let Y 1 (x,t),Y 2 (x,t) be two branches of 
the solution Y(x, t) in O, that is the functions Y\(x, t), Y 2 (x, t) are obtained 
by the analytic continuation of the solution of (|6.1|) along two different paths 
71,72 leading from x to O. These branches Y\(x, t), Y 2 (x, t) must satisfy 
the relation Y 2 {x, t) = Y\{x, £)M 7 (i), where M 7 (£) is the monodromy matrix 
corresponding to the loop 7 constructing from the paths 71 and 72. In gen- 
eral, this monodromy matrix M 7 (i) may depend on t. However, according 
to Theorem 17.121 the considered deformation is isomonodromic. In partic- 
ular, the matrix M 7 (£) does not depend on t, that is Y 2 (x, t) = Yi(x,£)M 7 , 
where M 7 does not depend on t and is invertible. (M 7 may depend on the 
choice of the branches Y\ and Y 2 .) Differentiating the last equality with 

QY 2 dY\ 

respect to t k , we obtain — — = — — M 7 From the last two equalities it fol- 
ate Ot k 

9Y\ dY 2 dY 

lows that — — YT = -r. — , that is the logarithmic derivative — — Y^ 1 

dt k 1 dt k 2 & _ dt k 

is a single- valued holomorphic function of x in C \ {t\, . . . , t n }. Using the 
representation (|7.2|) in a neighborhood of the singular point tj, we see that 
for j 7^ k 



?L({x,t))Y-\(x,t)) = dHa ^ t> H-j{x,t), x G D t . jP . 
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(For j ^ k, the second factor E a .(x — tj) in (|7.2|) does not depend on 
ifc.) Since the function H aj (x,t) is holomorphic with respect to (x,t) for 
x G D t . iP ,te D(t°,6), so is the function f£((x, t))Y -1 ((a:, t)) for j ^ fc. 
For j = A;, we have to take into account the factor E aj (x — tj): 

BY 

— (( x ,t))Y- 1 ((x,t)) = 

= H ak (x,t) {Jf k E ak {x - t k )E-*{x - t k )j H~l(x,t) 

+ 9Ha £' t) H£(x,t), xeD tk , p . 

dH a (x,t) 

The function ^ — - — H~Hx,t) is holomorphic for x € Df j:P . Taking into 

account (|8.2|) and (|8,5|) . we come to the conclusion: 

(8.6) ^(( Xl t))Y- 1 ((x,t)) = —^-+ 
at k x — tf. 



+ (a function holomorphic with respect to x), x £ -D^ 



fc,p- 
<9Y 

Thus, the function — — ((x,t))Y~ 1 ((x,t)), considered as a function of the 
dtk _ 

variable x, is holomorphic in the entire extended complex plane C, except 

the point t\~. Near tk this function behaves as it described in (|8.6|) . Therefore, 

^(( x ,t))Y-H(x,t)) = —\ + C k (t), 
dt k x -t k 

where Ck(t) does not depend on x. To determine C&(t), we use the normal- 
izing condition 15151 . Since Y(x ,t) = I, dY ^o,t) = o, i.e. - ® k ^ + 

dt k x-t k 
C k (t) = at x = x . Thus, <JHHJ> holds. Q.E.D. 

Theorem 8.2. Let akxk matrix functions Qj(t), t = (t%, . . . ,t n ), j = 
l,...,n, are holomorphic for t G T>, where T> C C™ is an open connected 
set, and satisfying the condition 

Qj(t) = 0, t£T>. 

l<j<n 

Assume that the deformation 

dx \lt<n X -^) 

is isoprincipal (in the wide sense) with respect to the distinguished point 
xq = oo at least for one pole loci t° £ T>. 

Then the matrix functions Qj(t),j = 1, . . . , n, satisfy the Schlesinger sys- 
tem TT1\) for all t£T>. 
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PROOF. Let Y(x,t) be the solution of the differential equation 1)6.1)1 satis- 
fying the normalizing condition Y(x,t) x=00 = I. According to Lemma 18,1) 
the equations (|4.1b|) are satisfied for x G C \ {t l5 . . .,t n },t G D(t°,<5). The 
equation (|4.1aj) is the same as the equation 1)6. According to Theorem 
14.11 the matrix functions Qj(t), t = (ti, . . . ,t n ), j = 1, ...,n satisfy the 
Schlesinger system p. 1)1 for t € D(t°, 5). Since matrix functions in both left 
and right hand sides of 1)1-1)1 are holomorphic in X> and coincide in D(i°, 5), 
they coincide everywhere in T>. Q.E.D. 

9. CONSTRUCTION OF THE ISOPRINCIPAL DEFORMATIONS 
AS A TOOL FOR SOLVING THE SCHLESINGER SYSTEM. 

Theorem 18.21 opens a way for constructing solutions of the Schlesinger 
system. The Cauchy problem for this system can be formulated as follows. 
Given t° G C™ and k x k matrices . . . , Q^, one needs to find matrix 
functions Qi(t), . . . , Q n (t) of the variable t satisfying the Schlesinger system 
i)Li|) and the initial condition 



enumerated by the pole loci t = (ti, . . . ,t n ), with holomorphic coefficients 
Qi(t), ... ,Q n (t) (a holomorphic deformation) such that this holomorphic 
deformation is isoprincipal with respect to xq = oo, then, according to 
Theorem 18.21 these coefficients will satisfy the Schlesinger system. In order 
for the initial conditions ()9.1)) to be satisfied, the family must contain the 
equation 



which corresponds to the parameter value t = t . 

If all the matrices Q®, . . . , Q„ are non-resonant (and hence, according to 
Theorems 17.121 and 17.51 all the matrices Qi(t), - - - , Q n (t) are non-resonant 
as well), this deformation is isoprincipal if, and only if, it is isomonodromic 
(Theorem 17.121 and Lemma 17. 7)1 . Therefore, if all the matrices Q?, • • • , Qn 
are non-resonant, the Cauchy problem with initial conditions ()9.1)) is reduced 
to the problem of constructing the isomonodromic deformation ()9.2)) which 
in particular contains the equation ()9.3)1 , An approach to constructing such 
isomonodromic deformation lies in using the Riemann-Hilbert problem. 

The Riemann-Hilbert problem can be formulated as follows. 
Given t = (ti,...,t n ) G C™, loops 71 (t), - - - , 7n(t) with the distinguished 
point xq that form the geometric configuration as it is plotted in Figure Q 



(9-1) Q j (t°) = Q° j , j = l,...,n. 



If we succeed in constructing a family of Fuchsian equations 



(9.2) 




(9.3) 
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and matrices Mi, . . . ,M n satisfying the condition M^m ^ir(n) = !• 
One needs to construct the Fuchsian differential equation i)R2)J with poles 
ti, . . . , t n whose monodromy matrices take the prescribed values 

(9.4) M 7 . (t) = M,-, j = l,...,n. 

An approach for solving the Riemann-Hilbert problem was developed by 
J. Plemelj in 1908, and, a little later but independently, by G.D. BirkhofT 
(see [Plemlj . |Plem2j . [Birkj ) . The first step of this approach requires to 
solve a factorization problem for a matrix function on a closed contour. The 
contour passes through points ti,...,t n , and the matrix function is piece- 
wise constant, constructed from the monodromy matrices Mi, . . . , M n . This 
factorization problem can be reduced to a singular linear matrix integral 
equation on the contour. On the second step the solution of the factorization 
problem undergoes certain additional transformations in order to obtain 
the solution of the original Riemann-Hilbert problem. The second step is 
not completely painless: there are reefs here. In fact, not long ago A. A. 
Bolibruch has discovered examples of data for which the Riemann-Hilbert 
problem has no solution (see |AnBoj . [Boll , Bol2 j, BokfJ). 

In order to apply the Riemann-Hilbert problem for the construction of 
the isomonodromic deformation, first we have to determine for the equation 
(|9.3|) (whose coefficients - matrices Q\ , . . . , - are the initial values in 
the Cauchy problem for the Schlesinger system) the monodromy matrices 
Mi, . . . , M°, corresponding to the loops 71 • • • , j n (t ). Then we have to 
solve the Riemann-Hilbert problem for the following data: an arbitrary point 
t = (ti, . . . , t n ) € C", the appropriate loops 71 (i), • • • , 7n(*)> and the matri- 
ces Mi, . . . , M°. As a result we shall obtain the matrices Qi(t), . . . , Q n (t). 
Repeating this procedure for every t, (and, therefore, solving the infinite 
family of the Riemann-Hilbert problems) we construct matrix functions 
Qi(t), . . . ,Q n (t). If we succeed in performing this construction in such a 
way that these matrix functions are holomorphic then we obtain the desired 
isomonodromic deformation. Finally, if the initial data are non-resonant, we 
obtain the solution of the Cauchy problem for the Schlesinger system. Thus 
solving the Cauchy problem for a non-linear system is reduced to solving a 
family of linear problems. 

However, solving the Schlesinger system in the non-resonant case is not 
the ultimate goal of the present manuscript. We would like to consider 
the case which is, in a certain sense, just the opposite. More precisely, we 
consider the initial values Q\, . . . , Q° such that fundamental solutions of the 
Fuchsian system (|9.3[) are rational matrix functions of x in general position. 
Of course, a rational function is uni-valued in the complex plane, hence 
the monodromy in this case is trivial: M\ = . . . = M n = I. In this case 
the spectrum of each matrix Q® is either {0, 1} or {0, —1} and rank(Q^) = 

^{tti, . . . , 7r n } is a permutation of the indices {1, . . . , n} which is determined by the geometric 
configuration of the paths a\, . . . , a n . For the configuration plotted in Figure [Tithe permutation 
is trivial: = k. 



48 



V.KATSNELSON AND D. VOLOK 



1. In particular, all the matrices Q®, . . . are resonant. Therefore, the 
considerations above, concerning the isomonodromic deformations, are not 
applicable in this case. Moreover, any family of Fuchsian systems with 
rational solutions is isomonodromic (the monodromies are trivial). Thus 
there are "too many" isomonodromic deformations here, and only one of 
them leads to the solution of the Schlesinger system. This is the isoprincipal 
deformation. 

Let us recall Theorem 15.161 and, in particular, representation 1)5. 28|) . If 
system (|5.1[) has a generic rational solution Y(x) then in 1)5.28(1 we have 
Aj = 0, and hence Y(x) is locally representable in a neighborhood of x = tj 

as 

(9.5) Y(x) = Hj(x)(x -tjfi, 
with Hj(x) holomorphic and invertible there, and 

(9.6) Z] = ±Zj. 

Furthermore, under condition ()9.6J) the factor (x — tj) Zj has the following 
simple form: 

( I-Zj + (x- tj)Zj, if Z] = Zj, 

( 9 - 7 ) - ^ = I + Zj — -*—Z it if Z* = -Zj. 

y x — tj J 

Thus the factorization (|9.5|) takes the form 
(9.8) Y(x) = Hj(x) (i - Zj + {x- tj)Zj) , if Z) = Z } , 

2 



(9.9) Y(x) = Hj(x)[l + Zj Zj), if Zj = -Zj. 

Factorizations of the form H9.8|) , ()9.9|) are traditional in the theory of rational 
matrix functions. The factors (|9,7|) are just the principal factors of the 
rational solution Y{x) at the singular point tj. If Y{x,t) is a family of 
rational solutions for a family of Fuchsian systems which gives an isoprincipal 
deformation then for every t = (ti, . . . ,t n ) at each singular point x = tj 
Y(x, t) must admit the factorization 



j > 



(9.11) Y{x,t) = Hj{x,t)[l + Zj - — ~—Zj), if Z) = -Zj, 



(9.10) Y(x, t) = Hj(x, t) [I - Zj + {x- tj)Zj ) , if Z) = Z 

1 

where Zj do not depend on t, and Hj(x,t) are holomorphic and invert- 
ible at x = tj. The realization theory for rational matrix functions allows to 
construct families of Fuchsian systems whose solutions Y (x, t) admit the fac- 
torizations (|9.10|) , (|9.11j) . Development of the realization theory appropriate 
for this goal originated with L.A. Sakhnovich [Sakhnj . This topic was fur- 
ther developed in |GKBMj . In [Katslj and |Kats2j this theory was adapted 
for applications concerning the Schlesinger system. The forthcoming second 
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part of the present manuscript will be dedicated to detailed construction of 
rational solutions of the Schlesinger system, based on the realization theory 
for rational matrix functions. 
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AVVEM'DTX 



A. THE REPRESENTATION OF A SOLUTION Of A FUCHSIAN 
EQUATION IN A NEIGHBORHOOD OF A SINGULAR 
POINT IN THE NON-RESONANT CASE. 

PROOF OF THE PROPOSITION O Choose a positive 6 so small that the 
closure of the polydisk D(t°, 5) is contained in the set T>, where the matrix 
functions Q p (i),p = 1, . . . , ra, are defined and holomorphic. We also impose 
the condition 5 < |pj(t°). So, if \tj — tj\ < 5, and \t p — t®\ < S,p ^ j, 
then \tj — t p \ > pj(t°). Since the matrix Qj(t°) is non-resonant, choosing, 
if necessary, a smaller 5, we can also assume that the matrix Qj(t) is non- 
resonant for all t € D(t°,(5). Moreover, if 5 > is small enough, then the 
differences of eigenvalues A p (t) — X q (t) of the matrix Qj(t) are separated 
from non-zero integers: 
(A.l) 

\X p (t)-X q (t)-m\ > e > Vp,q = l,...,k, Vme Z\{0}, Vt € D(t°,5), 

where e > does not depend on t € D(t°, S), m G Z \ {0}. 
Changing variable 

(A. 2) x->x + tj, 

we reduce the differential equation (|6.1|) to the form 

(A.3) ^(x,t) = /Q 2 (t) + \ 

ax \ x J 

where 

(A.4) Z(x,t)=Y(x + tj,t), *(ar,t)- ^ Qp(t) 



x -{- tj i E 

<p<r 
V+3 



l<p<n ' ^ ^ 



In view of the relation \tj — t p \ > pj(t°), the function <£(x,i) admits expan- 
sion of the form 

(A.5) $(x,t) = $r(t)x r , xeC, I 

0<r<oo 

where are fcxfc matrix functions, which are holomorphic with respect 

to t for t £ D(i°,<5) and admit the estimate 

(A.6) ||$r(*)|| <Ci-(p j (t°))- r , teD(t°,8), 0<r<oo, 

where Ci < co does not depend on t and r. 

We seek a normalized solution of (jA.3|) . corresponding to the singular 
point x = 0, in the form 

(A.7) Z(x) = ^(x,t)x Qj{t) , 
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where \I/(x,i) is power series in x: 

(A.8) V(x,t) = ^2 ^At)x r , with W (t)=I. 

0<r<oo 

Substituting the expressions (lA~7l) - <TP|) and (TOl) into (TOl) and using 
the differentiation rule 

(A.9) — x^W = ^a^O, 

ax x 

we obtain for r = 0, 1, 2, . . . the equalities 



(A.10) (r+lJ^r+i^+^r+iCtjQ^CtJ-QjCtJ^r+iCt) = *,(*)*tn(t), 

i+m=r 
Z>0,m>0 

relating the given sequence {3> r (*)}o<r<oo and the sequence {\? r (*)}o<r<cx>) 
which has to be found. The relations (jA.lOj) express the fact that (|A.7|) - 
(|A.8|) is a formal solution of the differential equation (|6.1|) . If the series (|A.8|) 
converge in a disk {x G C : \x\ < p}, with p > 0, then (|A.7|) . where ^(x, t) is 
the sum of these series, is an actual solution of (|6.1|) in the punctured disk 

{x S C : |x| < p} \ {0}. 

Since the constant term ^o(t) is given: ^o(t) = I: the equations (|A.10|) 
can be considered as a recursive system for successive determining of *S>i(t), 
^2 (*) 5 etc. The r-th equation of the system (|A.10|) is of the form 

(A.ll) XX-[Q,X]=Y, 

where [Q, X] = QX — XQ — the commutator of the matruces Q and X. In 

(HZEqJi, 

(A.12) A = r + 1, Q = Qj{t), X = * r+1 (t), Y= *I (*)*»«(*)■ 

i+m=r 
Z>0,m>0 

Now we interrupt the proof of Proposition 16. 31 to discuss the matrix equa- 
tion (|A~TT|) . 

Lemma A.l. Lei 971/% 6e i/ie set of all k x k matrices with complex entries. 
Given a matrix Q £ 9Jtk, 1st us associate with this matrix the operator 



(A.13) ad Q X = QX - XQ, ad Q :Wl k ^VJl k . 

Let Ai(Q), . . . , Afc(Q) be the set of all eigenvalues of the matrix Q (enumer- 
ated with multiplicities). Then the set of all eigenvalues (enumerated with 
multiplicities) of the operator adg is the set of all the differences {A P (Q) — 

^q(Q)}l<p,q<k- 

PROOF. Assume first that all eigenvalues of the matrix Q are pairwise dif- 
ferent. Let Up and v q be, respectively, eigenvector columns and eigenvector 
rows of the matrix Q: 



Qu p = XpUp, V q Q = XgVg, 1 < p, q < k. 
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Then the matrices X pq = u p ® v q are eigenvectors of the operator adg: 

a,d Q X pA = (X p - X q )X Pjq , 1 < p, q < k. 

Thus, the numbers X P (Q) — X q (Q) are eigenvalues of the operator adQ. Since 
the eigenvalues of the matrix Q are pairwise different, the vectors {u p }i< p <k 
and {v q }i< q <k form bases in the spaces of all /c-vector columns and fc-vector 
rows respectively. Since dimension of the space !J% is equal to k 2 , the set 
of matrices {Xp,g}i< Pi(? <fc is a basis of the space Wl^- Therefore, the set 
{X P (Q) — Xq(Q)}i< Piq <k is the set of all eigenvalues of the operator adQ 
(enumerated with multiplicities). This result can be extended to general 
matrices Q, i.e. to matrices whose eigenvalues are not necessarily pairwise 
different. Such an extension can be done using the approximation reasoning. 
Q.E.D. 



From the above result it follows that if A satisfies the condition A 7^ X P (Q) — 

1 < Pi 1 < fe, then the equation (Q, X, Y are k x k matrices) 
(A. 14) XX - ad Q X = Y 

is solvable with respect to X for every Y. However, for what follows one 
needs not only the solvability of this equation but also an estimate for its 
solution. 

Lemma A. 2. Let C be a K -dimensional vector space equipped with a norm 
\\.\\c, and let A be a linear operator acting in C. Assume that the eigenvalues 
A m (-4) of the operator A are separated from zero: 

(A.15) |A m (,4)|>£, l<m<K. 

Then the inverse operator A^ 1 admits the estimate from above: 

(A.16) 



\A-H <2 K \\A\\ K - 1 e- K , 



where \\A\\ = sup\\Ax\\c,\\A 1 | 
all x € C such that \\x\\c < 1. 



sup \\A 1 x\\c, and sup's are taken over 



PROOF. Let x(C) = det(£Z — A) be the characteristic polynomial of the 
operator A, 

x(C) = C K - ^iC^ 1 + • • • + (-I)*" V-iC + {-i) K <T K , 

where 07, 1 < / < K, denotes the l-th elementary symmetric polynomial of 
in Xi(A), ... , 

Xk(A). According to Hamilton-Cayley theorem, x(~A) = 0, hence 

0~K 0~K 

^From ()A.17|) it follows that 

K-2 , , °K-\ 



(A.17) A 



-1 



-1 



< 




\\A\\ K - l + 






0~K 







+ 



+ 



(*0 = 1) 
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It is clear that the fractions — are the elementary symmetric polynomials 

o K 

in the inverse eigenvalues - — -—- ..... - — ttt • From (IA.15I) we conclude that 

Ai(^l) X K (A) 



< C l K e l ~ K , where C l K are binomial coefficients: C l K 



l\{K-l)V 



Therefore, 



II^H < C Q K e- K \\A\\ K - 1 + C l K e l - K \\A\\ K - 2 + ■■■ + Cf-V 1 
Using the binomial formula, we rewrite the previous inequality as 

(A.18) \\A- l \\ < \\AW~ 1 ((1 + £ -1 |L4||) K - 1 

Since |A m (.4)| < \\A\\ for any eigenvalue X m (A), inequality (|A.15|) implies 
that 1 < e" 1 !!^ !!, and hence 1 + < 2e~ 1 ||„4||. Weakening the in- 

equality (|A.18|) . we come to the more rough, but more simple inequality 
<fA~Tfil) . Q.E.D. 

Lemma A. 3. Let Q be a k x k matrix, with eigenvalues Ai,...,Afc, enu- 
merated with multiplicities. Let X be a complex number, and let for some 
e > the inequalities 

(A.19) |A-(A P -A,)| >e, l<p,q<k, 

hold true. Then for every Y, the equation \A. 1J$ ( or, what is the same, the 
equation \A.11\) ) is solvable with respect to X, and the estimate 

(A.20) \\X\\ < 2 fe V fe2 (|A| + 2||Q||) fc2 - 1 ||Y||, 

holds, where \\.\\ is an arbitrary norm on *DTfc with the property: ||^4-B|| < 
11^1111^11 for every A,B € Wl k . 

PROOF. We use Lemma lA.2l for the operator A = XI — adg, acting on the 
space 9Jtfc. It is clear that ||^4|| < |A| + 2||Q||. According to Lemma TA.ll the 
set {A — (A p — A 9 )}i< p . 9 <fc is the set of all eigenvalues of A (enumerated with 
multiplicities). The inequality (jA.19|) provides an estimate of these eigen- 
values from below. Applying Lemma IA.1I to this operator A (and taking 
into account that dim9Jlfc = k 2 ), we obtain the estimate (|A.20|) . Q.E.D. 



The following lemma gives an estimate, which is more precise than (TA~2"0l) 
for large A. 

Lemma A. 4. Let Q be a k x k matrix, with eigenvalues Ai,...,A&, enu- 
merated with multiplicities. Let X be a complex number which is different 
from the numbers X p — X q , 1 < p, q < k, and let the estimates for \\Q\\ from 
above and for \X — (X p — X q )\ from below be given: 



(A.21) |A-(Ap-A 9 )| > e, l<P,q<k, 



\\Q\\ < H, (e > 0, n < oo). 
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Then for every Y £ 9Jik, the equation \A.lJj\) (or, what is the same, the 
equation \A.11\) ) is solvable with respect to X € 97tfc, and the estimate 

(A.22) 11X11 < C(£ '^ fe) IIFIL 

\ ) ii ii - 1 + | A | ii ii, 

holds, where J.J is an arbitrary norm on with the property: \\AB\\ < 
||-A||[|JB|| for every A,B£ 9Jtk, and C(e,fi,k) < oo is a constant depending 
only on e, /i and k. 

PROOF. If XX - ad Q X = Y, then |A|||X|| < ||ad Q X|| + ||Y|| < 2||Q||||X|| + 
||Y|| < 2fi\\X\\ + Thus, 

(A.23) < — J_||y||, for |A|>2/x. 

On the other hand, if we assume that |A| < 3/i and replace in the inequality 
HA.20|) |A| and ||Q|| by the larger values 3// and \i respectively, we come to 
the inequality 

(A.24) ||X|| < W^e'^fi^^WYW, for |A| < 3/i. 

Unifying the inequalities ()A.23|) and ()A.24|) . we come to the inequality (|A.22|) 
with a suitable constant C(e,fi,k) (which, of course, can be found explic- 
itly). Q.E.D. 

Now, after we have investigated the matrix equation (jA.llj) . we resume 

the proof of Proposition lb. 31 We apply Lemma IA. 31 to the equality (|A.10|) 

considered as an equation of the form (|A.11|) . with X, Y, Q,X defined in 

(lA~12l) . The numb er max ||Qj(t)|| serves as \x, the number e is taken from 
teD(t°) 

(|A.1|) . Substituting these X, A and the estimate 

(A.25) ||Y|| < Yl ll*r-m(*)||||*m(*)|| 

0<m<r 

for \\Y || into ()A.22|) . we obtain, for t € D(t°, S), r = 0, 1, 2, . . . the estimate 

(A.26) H* r+ l(t)|| < £ ||*r-m(*)||||* m (*)||. 

0<m<r 

Substituting the upper estimate (|A.6|) for ||$ r _ m (t)|| into (|A.26|) . we obtain: 



(A.27) 

\\* r+1 (t)\\p r+1 < E \^rn(t)\\P m , VteD(t°,5), r = 0,l,2,..., 

0<m<r 



r + 2 
where 

(A.28) C 2 = p-C 1 -C(e,n,k) 
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p = pj(t : ) is defined in (|fi.7j) . Ci and C{e,p,,k) are the values from ()A.6|) 
and ()A.22|) respectively. It is important to note that the value C2 < 00 does 
not depend on r and on t € D(i°). 

Lemma A. 5. Let {a r }o< r<00 be a sequence of non-negative numbers, sat- 
isfying the conditions 

(A.29) a r+1 < — ^— dm, r = 0,1,2,...; a = 1, 

0<m<r 

where d, < d < 00, does not depend on r. 

Then the sequence {a r }o< r<00 admits the upper estimate 

(A.30) a r < (r + l) d_1 , r = 0,1,2,.... 

PROOF. For r = 0, the assertion (|A.30|) is true: ao = 1. Assume that the 
inequalities (|A.30|) hold for r = 0, 1, . . . ,m. Then, according to (|A.30|) . the 
inequality 

(A.31) a m+1 < -L- (r + l) d 1 

m + 2 f-T 

U<r<m 

holds. If d > 1, then 

r+2 m+2 



(r + l)*" 1 < J t^dt, Yj ( r + ^ / td ~ ldt ^ 



[m + 2) 



r 1 j 0<r<m 



and the inequality (|A.30|) holds for r = m + 1. If < d < 1, then 

r+1 m+1 



(r + l)^ 1 < J t^dt, Y ( r + ^ / td ~ ldt ^ 



0<r<m 



(m + 2) c 
d 



and the inequality (|A.30|) holds for r = m + 1 as well. By induction, the 
inequality (|A.30|) holds for all r. Q.E.D. 



Applying Lemma fA.5l to QA.27|) . with a r = \\^f r (t)\\p r , we obtain the estimate 



(A.32) ||*r(*)|| < {r + l) d p~\ r = 0,1,2,... , Vt€D(t°,5), 

where p = Pj(t°) is defined in (|6.7|) . d = C2 from (|A.28|) . d < 00 does not 
depend on r and on t € D(i°, <5). 

^From the estimate ()A.32|) it follows, that the power series (|A.8j) converge 
locally uniformly in the open circle {x : \x\ < p}, and the convergence is 
uniform with respect to the parameter t E D(t°,<5). Therefore, the sum 
fy(x,t) of these power series is a function, holomorphic with respect to x,t 
for {x : \x\ < pj(t )}, t € D(t°,<5). Let us prove that the function \P 



56 V.KATSNELSON AND D. VOLOK 

is holomorphic for these well. Since *f?o(t) = I, we can consider the 

formal power series 

(A.33) ¥ m {x,t)= ^r nv (*K, 

0<r<oo 



which are inverse to the series (|A.8|) . i.e. 
(A.34) 

( E Mt)x r ) ■ { E *r(t)z r ] 

V 0<r<oo J V 0<r<oo / 

= ( E ^ DV (t)x r ) • ( E *r{t)x r ) =1. 

yO<r<oo J yO<r<oo J 

It is clear that 

(A.35) ^ nv (t) = /. 

The expressions (|A.3|) . (|A.7|1 . ()A.8|) and (|A.34|) imply that the product 

(A.36) Z inv (x, t) = x- Q i^¥ nv (x, t) 

is a formal solution of the differential equation 



(A.37) — — = -Z mv + $(x,t) , 

ax y x / 

As before (see (|A.10j) ). substituting the expression (lA~36l) into (lA~37|) . we 

obtain the infinite recursive system 

(A.38) 

(r+l)*S 1 (t) + *^ 1 (t)Q i (t)-O i (t)*^ 1 (t) = - 2 *£W(t). 

i+m=r 
£>0,m>0 

Like the r-th equation of the system ()A.10|) , the r-th equation of the system 
(|A~38)) is of the form (|A~TT|) . with the same A and Q as i n l|A.10|) : A = 
r + 1, Q = Qj(t). The only difference is that the right-hand side of (|A.11|) 

is Y = E**(*)*m(<) for (|A~TU1) . and Y = - E *m V (*)$/(*) f o r $EM- 
Because A and Q for both matrix equations, corresponding to the systems 
(|A.10|) and ()A.38|) . are the same, the estimates (jA,22|) for both equations 
hold with the same constant C(e,fi,k) (Of course, e,fi,k are the same). 
Instead of the estimate (|A.25j) . the estimate 

(A.39) ||Y||< ll*r-m(t)||||«^(t)|| 

0<m<r 



Of course, in view of the condition 'l'o(t) = I? the function 1 (x,t) is holomorphic for x 
which are small enough. However, it is very easy to prove the holomorphy of the function i S/~ 1 in 
the whole disc {x : \x\ < p}. Such a proof can be done "for free", using the previously established 
estimates. 
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for the right-hand side of (|A.11|) should be used now. Finally, we come to 
the recursive system of inequalities 

(A.40) < C{£ ^ k) £ \\^ m (t)\\\\¥^(t)\\, VteD(tV)- 

0<m<r 

with respect to || 1 J , ^ v (t)||. However, the latter system of inequalities coin- 
cides with the system (|A.26|) of inequalities with respect to ||vl/ m (t)||. Since 
||^ nv (t)|| = p (t )|| = || J||, th e estimates for ||tf m (t)|| and ||^ v (t)||, which 
follow from (|A.26|) and (|A.40|) . coincide. Thus, the estimates 

(A.41) ||*r(t)|| <(r + l)V r , r = 0,l,2,..., Vt€D(t°,5), 

hold with the same p and d that in ()A.32|) . Therefore, the series (jA.33|) 
converge locally uniformly in the disk {x : |x| < Pj(t)} and represent the 
matrix function ^ mv (x,t) which is holomorphic with respect to x,t for x : 
\x\ < pj(t°), t G D(t°, 5). The identities (TOD) for formal power series imply 
the identities 
(A.42) 

y(x,t)^ mv {x,t) = ¥ nv (x,t)V(x,t) = I V x : \x\ < pj(t°),t g D(t°,5) 

for the appropriate holomorphic matrix functions. Thus, the values of the 
matrix function ^(x,t) for x : \x\ < pj(t°), t S D(t°,<5) are invertible, and 
*- x (x,t) = ¥ nY (x,t). 

Performing the change of variables x — > x — tj, inverse to the change 
(|A.2f) . we obtain the representation (|6.3j) for the local solution Yj(x,t), with 
Hj(x,t)) = ^(x — tj,t). The properties of the function Hj(x,t), stated 
in the formulation of Proposition 16.31 express the properties of x , is) as 
established above. This completes the proof of Proposition 16.31 Q.E.D. 
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